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Figure  34  Po('*o’  Probability  Density 

Function  for  Successive  Zero  Crossing 
Intervals  of  a  Random  Process  Consisting 
of  Sine  Wave  of  Frequency  ^/2  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  WQ(f).  The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 
Theoretical  Approximations  are  Compared 
with  Experimental  Results. 

Figure  35  P^Cu  ,  a)  is  the  Probability  Density  Function 
for  the  Sum  of  Two  Successive  Zero- 
Crossing  Intervals  of  a  Random  Process 
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fo/2  Plus  Gaussian  Noise  Having  the 
Power  Spectral  Density  W  (f).  The 
Ratio  of  the  Average  Sine  Wave  Power 
to  the  Average  Noise  Power  is  Denoted 
by  the  Parameter  a. 

Figure  36  Pa(''»  ^^o)  the  Probability  that  a  Given 

Interval  Uq  Contains  Exactly  n  Zero- 
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fo/2  Plus  Gaussian  Noise  Having  the 
Power  Spectral  Density  Wo(f).  The 
Ratio  of  the  Average  Sine  Wave  Power 
to  the  Average  Noise  Power  is  Denoted 
by  the  Parameter  a. 

Figure  37  A  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
Wo(f)  Defines  Zj^{n,  u^),  the  Conditional 
Probability  That  the  nth  Zero-Crossing 
Point  From  a  Given  Zero-Crossing 
Point  in  duo  Occurs  After  the  time  duo+Ug. 
The  Ratio  of  the  Average  Sine  Wave 
Power  to  the  Average  Noise  Power  is 
Denoted  by  the  Parameter  a. 
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Figure  41 


Probability  Density  Function 
for  Successive  Zero-Crossing  Intervals 
of  a  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  i^/z  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
W'o(f)‘  The  Ratio  of  the  Average  Sine 
Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 
Theoretical  Approximations  are  Compared 
With  Experimental  Results. 

Pl(uQ,a)  is  the  Probability  Density 
Function  for  the  Sum  of  Two  Successive 
Zero-Crossing  Intervals  of  a  Random 
Process  Consisting  of  a  Sine  Wave  of 
Frequency  i^/l  Plus  Gaussian  Noise 
Having  the  Power  Spectral  Density 
WQ{f).  The  Ratio  of  the  Average  Sine 
Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 

P2^(n,  Uq)  is  the  Probability  that  a 
Given  Interval  Uq  Contains  Exactly 
n  Zero-Crossing  Points  of  a  Random 
Process  Consisting  of  a  Sine  Wave 
of  Frequency  fo/2  Plus  Gaussian 
Noise  Having  the  Power  Spectral 
Density  Wo(f).  The  Ratio  of  the 
Average  Sine  Wave  Pew  er  to  the  Average 
Noise  Power  is  Denoted  by  the 
Parameter  a. 

A  Random  Process  Consisting  of  a 
Sine  Wave  of  Frequency  Iq/Z  Plus 
Gaussian  Noise  Having  the  Power 
Spectral  Density  WQ(f)  Defines 
Zj^(n,  Uq).  the  Conditional  Probability 
that  the  nth  Zero-Crossing  Point 
From  a  Given  Zero- Crossing  Point 
in  duo  Occurs  After  the  Time  dUg+UQ. 

The  Ratio  of  the  Average  Sine  Wave 
Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a. 
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Figure  43 
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Figure  45 


Po(Uy,  a)  ie  the  Probainlity  Density  Function 
for  Succ  essive  Zero-Crossing  Intervals 
of  a  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  Iq/z  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
Wo(f).  The  Ratio  of  the  Average  Sine 
Wave  Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a.  Theoretical 
Approximations  are  Compared  with 
Experimental  Results. 

Plfuj^.a)  is  the  Probability  Density  Function 
for  the  Sum  of  Two  Successive  Zero- 
Crossing  Intervals  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency 
fjj/?.  Plus  Gaussian  Noise  Having  the 
Power  Spectral  Density  Wo{f).  The 
Ratio  of  the  Average  Sine  Wave  Power  to 
The  Average  Noise  Power  is  Denoted 
by  the  Parameter  a. 

Pj^(n,  u^)  is  the  Probability  that  a  Given 
Interval  u^  Contains  Exactly  n  Zero- 
Crossing  Points  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency 
fg/Z  Plus  Gaussian  Noi sc  Having 
the  Power  Specrral  Density  Wg(f). 

The  Ratio  of  the  Avc'rago  Sine  Wave 
Power  to  the  Average  Noise  Power  is 
Denoted  by  thc'  Parameter  a. 

A  Random  Process  Consisting  of  a 
Sine  Wave  of  Frequency  fg/Z  Plus 
Gaussian  Noise  Having  the  Power 
Spectral  Density  WQ(f)  Defines 
Z  (n,  Uq),  the  Conditional  Probability 
That  the  nth  Zero -Crossing  Point 
from  a  Given  Zero-  Crossing  Point 
in  dUg  Occurs  After  the  Time  dUg+Ug. 

The  Ratio  of  the  Average  Sine 
Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the 
Parameter  a. 
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Figure  49 


Po(uo,  a)  ie  the  Probability  Density 
Function  for  Successive  Zero-Crossing 
Intervals  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  Wo(f).  The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a.  Theoretical 
Approximations  are  Compared  with 
Experimental  Results. 

Probability  Density 

Function  for  the  Sum  of  Two  Successive 
Zero-Crossing  Intervals  of  a  Random 
Process  Consisting  of  a  Sine  Wave  of 
Frequency  i^/Z  Plus  Caussian  Noise 
Having  the  Power  Spectral  Density 
Wo(f)-  The  Ratio  of  the  Average  Sine 
Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 

p^(n,  Uq)  is  the  Probability  that  a  Given 
Interval  u^  Contains  Exactly  n  Zero- 
Crossing  Points  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency 
fo/2  Plus  Gaussian  Noise  Having  the 
Power  Spectral  Density  Wo(f).  The 
Ratio  of  the  Average  Sine  Wave  Power 
Lo  the  Average  Noise  Power  is  Denoted 
by  the  Parameter  a. 

A  Random  Process  Consisting  of  a 
Sine  Wave  of  Frequency  i^/Z  Plus 
Gaussian  Noise  Having  the  Power 
Spectral  Density  WQ(f)  Defines  Zj^(n,  u^), 
the  Conditional  Probability  that  the  nth 
Zero-Crossing  Point  from  a  Given  Zero- 
Crossing  Point  in  duo  Occurs  After  the  time 
dUg+UQ.  The  Ratio  of  the  Average  Sine 
Wave  Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a. 
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Figure  5  3 


is  the  Probability  Density  Function 
for  Successive  Intervals  Defined  by 
Adjacent  Stationary  Points  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency  i^/l 
Plus  Gaussian  Noise  Having  the  Power  Spectral 
Density  WQ(f).  The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 
Theoretical  Approximations  are  Compared 
with  Experimental  Results. 

M|(u^,  a)  is  the  Probability  Density  Function 
for  the  Sum  of  Two  Successive  Intervals 
Defined  by  Adjacent  Stationary  Points  of  a 
Random  Process  Consisting  of  a  Sine  Wave  of 
Frequency  i^/l  Plus  Gaussian  Noise  Having 
the  Power  Spectral  Density  WQ(f).  The 
Ratio  of  the  Average  Sine  Wave  Power  to  the 
Average  Noise  Power  is  Denoted  by  the 
Parameter  a. 


m  (n,  u  )  is  the  Probability  that  a  Given 
Interval  Uq  Contains  Exactly  n  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  WQ(f).  The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 


A  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  f^/Z  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
Defines  Sa(n,  u^),  the  Conditional 
Probability  that  the  nth  Stationary  Point 
from  a  Given  Stationary  Point  in  dUo 
occurs  after  the  Time  du^+Uo.  The  Ratio 
of  the  Average  Sine  Wave  Power  to  the 
Average  Noise  Power  is  Denoted  by 
the  Parameter  a. 
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M  (Uq,  a)  is  the  Probability  Density  Function 
for  Successive  Intervals  Defined  by  Adjacent 
Stationary  Points  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency  Iq/Z 
Plus  Gaussian  Noise  Having  the  Power 
Spectral  Density  Wo(f).  The  Ratio  of 
The  Average  Sine  Wave  Power  to  the 
Average  Noise  Power  is  Denoted  by  the 
Parameter  a.  Theoretical  Approximations 
are  Compared  with  Experimental  Results. 

M|(uo,  a)  is  the  Probability  Density  Function 
for  the  Sum  of  Two  Successive  Intervals 
Defined  by  Adjacent  Stationary  Points 
of  a  Random  Process  Consisting  of 
A  Sine  Wave  of  Frequency  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  Wo(f).  The  Ratio  of  the 
Average  Sine  Wave  Power  to  the  Average 
Noise  Power  is  Denoted  by  the  Parameter 
a. 

ma(n,  uo)  is  the  Probability  That  a  Given 
Interval  Ug  Contains  Exactly  n  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  fg/Z  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  ''H The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 

A  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  fg/Z  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
Wg(f)  Defines  Sj^(n,  Ug),  the  Conditional 
Probability  that  the  nth  Stationary 
Point  from  a  Given  Stationary  Point 
in  du  Occurs  after  the  Time  du^+u^. 

The  Ratio  of  the  Average  Sine  Wave 
Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a. 


Figure  58  Mq(uo,  a)  is  the  Probability  Density  Function 
for  Successive  Intervals  Defined  by 
Adjacent  Stationary  Points  of  a  Random 
Process  Consisting  of  a  Sine  Wave  of 
Frequency  fo/Z  Plus  Gaussian  Noise 
Having  the  Power  Spectral  Density  Wo(f)  • 
The  Ratio  of  the  Average  Sine  Wave 
Power  to  the  Average  Noise  Power  is 
Denoted  by  the  Parameter  a.  Theoretical 
Approximations  are  Compared  with 
Experimental  Results. 

Figure  59  Mi(Uq,  a)  is  the  Probability  Density 

Function  for  the  Sum  of  Two  Successive 
Intervals  Defined  by  Adjacent  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  f  /z  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  W  ff).  The  ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parauneter  a. 

Figure  60  ma(n  ,  u^)  is  the  Probability  that  a  Given 
Interval  Uq  Contains  Exactly  n  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  f^/Z  Plus 
Gaussian  Noise  Having  the  Power 
Spectral  Density  Wo(f).  The  Ratio  of 
the  Average  Sine  Wave  Power  to  the 
Average  Noise  Power  is  Denoted  by 
The  Parameter  a. 

Figure  61  A  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  f^/Z  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
Wo(f)  Defines  Sg(n,  Uq),  the  Conditional 
Probability  that  the  nth  Stationary 
Point  from  a  Given  Stationary  Point 
in  duo  Occurs  adfter  the  Time  duo+u^  . 

The  Ratio  of  the  Average  Sine  Wave 
Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a. 


Figure  62 


Figure  63 


Figure  64 


Figvire  65 


a)  is  the  Probability  Density  Function 
for  Successive  Intervals  Defined  by  Adjacent 
Stationary  Points  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency  fo/2 
Plus  Gaussian  Noise  Having  the  Power  Spectral 
Density  Wo(f).  The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a.  Theoretical 
Approximations  are  Compared  with 
Experimental  Results. 

Mj^(uo.  a)  is  the  Probability  Density  Function 
for  the  Sum  of  Two  Successive  Intervals 
Defined  by  Adjacent  Stationary  Points  of 
A  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  f^/Z  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
V/q{{).  The  ratio  of  the  Average  Sine 
Wave  Power  to  the  Average  Noise  Power 
is  Denoted  by  the  Parameter  a. 

)  is  the  Probability  that  a  Given 
Uq  Contains  Exactly  n  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wav'o  of  Frequency  f^/Z  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  \VQ(f),  The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 


ma(n,  u 
Interval 


A  Random  Process  Consisting  of  a  Sine 
Wave  of  Frequency  f^/Z  Plus  Gaussian 
Noise  Having  the  Power  Spectral  Density 
Wo(f)  Defines  Sjj(n,  u^).  The  Conditional 
Probability  that  the  nth  Stationary  Point 
from  a  Given  Stationary  Point  in  du^ 
Occurs  After  the  Time  du„+u  .  The 
Ratio  of  the  Average  Sine  Wave  Power 
to  the  Average  Noise  Power  is  Denoted 
by  the  Parameter  a. 


Figure  66  Each  Intensity  Pattern  Represents  an 
Approximation  to  the  Joint  Probability 
Density  of  the  ith  and  (i+n)th  Zero- 
Crossing  Intervals  of  a  Gaussian  Proces 
Having  Power  Spectral  Density 

W3(f). 

Figure  67  The  Truncated  Gaussian  Probability 
Density,  Pq(t),  Represents  the  Most 
Random  Distribution  of  Successive 
Zero-Crossing  Intervals  Having  a 
Given  Mean  and  a  Given  Standard 
Deviation. 
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ABSTRACT 


This  report  describes  a  digital  system  for  studying  the  zero¬ 
crossing  intervals  of  random  procesctes.  Probabilities  and  probability 
densities  defined  by  the  zero- crossing  points  of  various  Gaussian 
processes  are  presented.  Probabilities  and  probability  densities  de¬ 
fined  by  the  zero-crossing  points  of  a  random  process  consisting  of 
a  sine  wave  plus  a  Gaussian  process  are  also  presented.  Finally, 
probabilities  and  probability  densities  defined  by  the  stationary  points 
of  a  random  process  consisting  of  a  sine  wave  plus  a  Gaussian  pro¬ 
cess  are  presented.  At  present  none  of  the  probabilities  or  probability 
densities  can  be  derived  explicitly  by  analytical  methods.  The  stan¬ 
dard  deviations  associated  with  the  probability  densities  are  also 
presented.  In  the  case  of  the  Gaussian  processes  the  correlation  co¬ 
efficients  for  two  successive  intervals  are  presented.  The  first  mo- 
mentsassociated  with  the  probability  densities  are  compared  with  the 
exact  theoretical  values.  All  the  other  experimental  results  are  com¬ 
pared  with  theoretical  approximations.  The  statistical  dependence 
between  the  ith  zero-crossing  interval  and  the  (i+n)th  zero-crossing 
interval  is  investigated. 


I. 


INTRODUCTION 


In  many  branches  of  science  and  technology  one  encounters  the 
basic  density  the  probability  density  of  the  interval  t  between 

the  mth  and  (m+n+l)th  zero-crossing  points  of  a  random  process. 

This  density  occurs  for  example  in  the  fields  of  statistical  commun¬ 
ication  theory,  oceanography,  statistical  mechanics,  biophysics, 
and  control  engineering.  When  n  =  0  the  interval  t  is  celled  a  zero- 
crossing  interval.  A  striking  experimental  fact  is  that  the  successive 
zero-crossing  intervals  of  an  audio  process  contain  a  great  deal  of 
the  intelligibility.  For  the  most  common  random  process,  a 
Gaussian  process  with  arbitrary  power  spectrum,  Pjj{t)  cannot  be 
derived  explicitly  by  analytical  methods.  For  a  Gaussian  process 
having  a  certain  class  of  power  spectra,  S.  O.  Rice  (1)  derived  two 
functions  which  can  be  used  to  approximate  TPq{t)  and  P^(t)  for 
small  T.  More  recently  J.  A.  McFadden  (Z)  and  M.  S.  Longuet- 
Higgins  (8)  derived  approximations  for  Pq{t)  which  compare 
favorably  with  experimental  results  for  a  larger  class  of  Gaussian 
processes  having  arbitrary  power  spectra. 

Some  other  important  theoretical  work  concerned  with 
has  been  reported  by  D.  Middleton  (3),  A.  Kohlenberg  (4), 

Kuznetsov,  Stratonovich,  and  Tikhonov  (5),  V.  I.  Tikhonov  and  I.  N. 
Amiantov  (6),  M.  S.  Longuet -Higgins  (7),  C.  W.  Helstrom  (9),H.  Steinberg, 
P.  M.  Schultheiss,  C.  A.  Wogrin,  andF.  Zweig  (10),  D.  S.  Palmer  (11), 

H.  Debart  (12),  D.  Slepian  (13),  W.  M.  Brown  (14),  J.  S.  Bendat 


-20- 


(15).  M.  Kac  (16),  I.  Miller  and  J.  E.  Freund  (17),  H.  Zuhrt  (29), 

S.  Ehrenfeld  (30)  et  al. ,  and  L.  I.  Bialyi  (31).  Fewer  analytical 
results  concerning  Pq(t),  P^(t)  exist  for  such  a  common  non- 
Gaussian  random  process  as  a  sine  wave  signal  plus  a  Gaussian 
process.  In  short,  one  must  conclude  that  mathematical 
difficulties  seriously  limit  the  understanding  of  this  statistical 
phenomenon.  Under  such  circumstances  one  is  naturally  motivated 
to  turn  to  experiment  for  further  insight. 

Most  of  the  experimental  work  dealing  with  this  problem  has 
been  reported  by  G.  M.  White  (18),  Favreau,  Low,  and  Pfeffer  (19), 
C.  R.  Gates  (20),  Kjell  Blotekjaer  (21),  and  A.  I.  Velichkin  and 
V.  D.  Ponomareva  (22).  In  order  to  explore  the  problem  further  an 
experimental  system  has  been  designed  at  this  Laboratory.  The 
system  differs  considerably  from  previous  experimental  systems, 
and  the  measurement  is  digital  rather  than  analog.  This  approach 
permits  the  system  to  work  in  conjunction  with  a  digital  computer 
to  help  analyze  the  recorded  data.  Also,  the  inherent  stability  of  a 
digital  system  allows  measurements  to  be  made  over  periods  of 
hours  or  days  if  required. 

This  report  presents  a  description  of  the  experimental 
system  and  the  following  results  concerning  the  zero-crossing 
intervals  of  certain  Gaussian  processes: 
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(1)  Pq(t),  the  probability  density  function  for  successive 
zero-crossing  intervals. 

(2)  P|(t),  the  probability  density  function  for  the  sum  of 
two  successive  zero-crossing  intervals. 

(3)  p(0,t),  the  probability  that  a  given  interval  t  contains 
exactly  zero  zero-crossings. 

(4)  p(1,t),  the  probability  that  a  given  interval  t  contains 
exactly  one  zero-crossing. 

(5)  the  conditional  probability  that  the  first 
zero-crossing  from  a  given  zero-crossing  in  dr 
occurs  after  the  time  dr  +  t. 

(6)  Z(2,t),  the  conditional  probability  that  the  second 
zero-crossing  from  a  given  zero-crossing  in  dr  occurs 
after  the  time  dr  +  t. 

(7)  The  expectations  and  standard  deviations  associated 
with  Pq(t)  and  Pj(t). 

(8)  The  correlation  coefficient  for  two  successive  zero¬ 
crossing  intervals. 

This  report  also  presents  the  following  results  concerning 
the  zero-crossing  intervals  of  a  random  process  consisting  of  a 
sine  wave  signal  of  frequency  plus  Gaussian  noise  such  that  the 
signal -to -noise  power  ratio  equals  a: 
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(1)  P^(Tta),  the  probability  density  function  for  successive 
zero-crossing  Intervals. 

(2)  P^(T,a),  the  probability  density  function  for  the  sum  of 
two  successive  zero-crossing  intervals. 

(3)  Pj^(0,t),  the  probability  that  a  given  interval  t  contains 
exactly  zero  zero-crossings. 

(4)  Pj^(f.T),  the  probability  that  a  given  interval  t  contains 
exactly  one  zero-crossing. 

(5)  Z^(i,T).  the  conditional  probability  that  the  first 
zero-crossing  from  a  given  zero-crossing  in  dr  occurs 
after  the  time  dr  +  t. 

(6)  Z^(2,t),  the  conditional  probability  that  the  second 
zero-crossing  from  a  given  zero-crossing  in  dr 
occurs  after  the  time  dr  +  t. 

(7)  The  expectations  and  standard  deviations  associated 
with  Pq(t,  a)  and  P^tr,  a). 

Finally,  this  report  presents  the  following  results  concerning 

the  intervale  defined  by  the  mathematical  stationary  points  of  a 

f 

random  process  consisting  of  a  sine  wave  signal  of  frequency 
plus  Gaussian  noise  such  that  the  signal-to-noiee  power  ratio 
equals  a: 

(f)  the  probability  density  function  for  successive 

intervals  defined  by  adjacent  stationary  points. 


-23- 


(2)  M|(T>a),  the  probability  density  function  for  the  sum  of 
two  successive  intervals  defined  by  adjacent  stationary 
points. 

(3)  m  (0,  t).  the  probability  that  a  given  interval  t  contains 

ft 

exactly  zero  stationary  points. 

(4)  m^(l.r)>  the  probability  that  a  given  interval  t  contains 
exactly  one  stationary  point. 

(5)  S^(I.t).  the  conditional  probability  that  the  first 
stationary  point  from  a  given  stationary  point  in  dr 
occurs  after  the  time  dr  4  r. 

(6)  S^(2,r),  the  conditional  probability  that  the  second 
stationary  point  from  a  given  stationary  point  in  dr 
occurs  after  the  time  dr  4  r. 

(7)  The  expectations  and  standard  deviations  associated  with 
M^(T.a)  and  M^(T,a). 

The  above  probabilities  and  probability  densities  were  selected 
for  investigation  not  only  because  they  are  of  theoretical  interest 
but  also  because  they  are  of  practical  interest,  since  they  are 
observables  which  are  comparatively  simple  to  measure. 

The  system  measures  the  basic  densities  P^(t).  P|(t). 

Po(Tta),  P|(T,a),  a)  and  M|(t,  a).  The  remaining  probabilities 

are  deduced  by  applying  the  theory  of  point-processes.  The  system 
can  measure  these  densities  at  levels  other  than  the  sero  level. 


II. 


THE  MEASURING  SYSTEM 


A,  METHOD  OF  MEASUREMENT 

Instead  of  measuring  successive  zero-crossing  intervals 
the  system  measures  samples  of  zero-crossing  intervals  which  are 
sufficiently  independent.  However,  the  random  processes  under 
investigation  are  metrically  transitive  (2  3)  or  ergodic,  and  hence 
the  statistics  of  these  samples  are  equal  with  probability  unity  to 
the  statistics  of  successive  zero-crossing  intervals. 

The  method  of  measuring  independent  samples  of  zero¬ 
crossing  Intervals  is  illustrated  in  Figure  1.  The  method  consists 
of  opening  and  closing  an  electronic  gate  at  the  proper  times 
following  each  initiate  command.  The  gate  controls  the  number  of 
clock  pulses  counted  by  an  electronic  counter.  The  resulting  count 
in  the  electronic  counter  is  then  a  measure  of  the  sampled  zero- 
crossing  interval.  The  count  is  printed  on  tape  in  order  to  monitor 
each  seimpled  time  interval  and  to  check  the  average  zero-crossing 
interval.  The  count  is  also  piinched  on  paper  tape  which  feeds  into 
an  IBM-7090  digital  computer.  The  circuit  details  have  been 
reported  in  the  literature  (24). 

Selecting  a  zero-crossing  interval  that  corresponds  in 
time  with  the  instant  of  an  initiate  command  tends  to  favor  the  long 
Intervals.  Since  in  general  successive  intervals  are  statistically 
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Initlolt  Commondt 


•—  Approiimotaly  I  tdc 


Figure  1 


Digital  Measurement  of  Independent 
Samples  of  Zero-Crossing  Intervals 
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dependent,  this  bias  may  also  bias  the  selection  of  an  interval  that 
occurs  after  an  initiate  command.  Accordingly,  the  digital  system 
skips,  with  equal  probability,  one  or  two  intervals  immediately 
following  the  initiate  command  and  then  measures  the  next  interval. 
By  skipping  these  intervals  the  bias  associated  with  selecting 
an  interval  after  an  initiate  conunand  proves  to  be  negligible. 

The  significance  of  this  bias  has  been  discussed  recently  by 
McFadden  (32).  The  system  has  provisions  for  skipping  an 
arbitrary  number  of  intervals. 

B.  SAMPLE  SIZE 

Consider  a  histogram  approximation  to  a  probability 
density  Pq(u)  as  shown  in  Figure  2.  After  measuring  n  mutually 
independent  samples  of  time  intervals  let  the  number  of  time 
intervals  in  the  range  u,  u  i-  6u  be  denoted  by  and  let 


u+6  u 


P^(u)  du 


Let  the  random  variable  Xj^  equal  one  if  the  result  of  the  kth 
measurement  is  in  the  range  u,  u  -f  6u  and  zero  otherwise.  Also, 
let  all  the  have  a  common  probability  density  with  mean  p  and 
variance  a^.  Then 


(2) 
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and  by  using  the  Tchebycheff  inequality 

Pr"^  |Sn-np|  ^  t  s  for  all  t  >  0  (3) 

where  Pr{E)  denotes  the  probability  of  the  event  E.  For  t  =  €n 


S 

Accordingly,  as  n  -►  oo,  -jp  -♦  p  with  probability  1.  This  is 
Bernoulli's  (i7i3)  celebrated  weak  law  of  large  numbers  which 
played  a  central  role  in  the  notion  of  probability.  Cantelli's  (1917) 
celebrated  strong  law  of  large  niunbers  not  only  Implies  the  weak 
law  of  large  numbers  but  also  asserts  that  with  probability  one  only 
finitely  many  of  the  events 
S 

1^  -  p|  >  €  occur  as  n  -*  oo. 

In  short,  the  weak  law  asserts  convergence  in  probability;  whereas 

the  strong  law  asserts  convergence  with  probability  one. 

In  experimental  work  n  is  limited  by  practical  considerations. 

Hence,  one  needs  a  quantitative  measure  of  the  manner  in  which 
S 

approaches  p  with  increasing  n.  In  deriving  such  a  measure 
notice  that  we  are  dealing  with  repeated  independent  trials  having 
only  two  possible  outcomes.  Accordingly,  we  will  adopt  a  Bernoulli 
model  and  consider  each  trial  as  a  Bernoulli  trial.  Using  the 
DeMoivre  limit  Theorem  (17i8)  and  letting  q  =  1-p  we  have  that  as 
n  -♦  00, 
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Pr{-J5jl^s  Ji  “PS  -x)=2  ®(x)-l 


(5) 


where 


*(x)  =  -iz.  J 

Zjt 


-y^/2 , 
e  '  dy 


-oo 


Let  US  require  that  -—  approach  p  to  within  a  given  arbitrary  amount 


n 


during  95%  of  our  atatlatical  experience.  That  la, let  our  confidence 


level  be  0.  95.  Then  x  =  1. 96, and  aa  ahown  in  Figure  2  la  given  by 


P  Au  =  I  -11  -  p  I  s  1 . 96 

o  '  n  '  N  n 


and 


AP  s  flf 


(6) 

(7) 


Histogram  Approximation  for  P^(u). 


Figure  2 
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Thla  ii  the  basic  result  concerning  the  sample  size.  It  was  used  to 
compute  the  95%  confidence  limits  presented  with  the  probability 
density  functions.  Using  this  basic  result  we  choose  n  to  be 
approximately  40,  000. 
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m.  THE  SYNTHESIS  SYSTEM 

A.  SHAPING  THE  POWER  SPECTRUM 

The  method  used  to  synthesize  a  Gaussian  random 
process  with  prescribed  power  spectral  density  is  illustrated  in 
Figure  3.  "White"  Gaussian  noise  is  applied  to  a  linear  network 
having  a  system  function  H(J2irf).  The  output  of  such  a  network 
is  a  Gaussian  random  process  having  a  power  spectral  density 
proportional  to  |HU2irf)|^.  Operational  amplifiers  with 
appropriate  feedback  networks  were  used  to  synthesize  H(jZTrf). 


Figure  3  Synthesis  of  Gaussian  Noise  with 
Prescribed  Power  Spectra. 


In  order  to  investigate  the  zero-crossing  intervals  of 
Gaussian  processes  having  certain  power  spectral  densities  of 
interest,  the  digital  system  requires  a  "white"  Gaussian  noise 
source  having  a  low  frequency  cut-off  of  approximately  1  cps  and 
a  high  frequency  cut-off  of  at  least  8  Kcps.  The  theory  and  design 
of  such  a  noise  source  has  been  reported  in  the  literature  (25). 

The  theory  is  described  below. 
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B.  THE  "WHITE"  NOISE  SOURCE 

Consider  a  train  of  periodic  pulses  u^^(t-nT^)  with  fixed 
duration  t  and  with  amplitudes  determined  by  the  corresponding 
amplitudes  of  an  ergodic  random  process.  This  random  pulse 
train  may  be  the  result  of  periodically  sampling  the  amplitude  of 
a  noise  voltage  y(t)  as  is  shown  in  Figure  4.  Generalizing 
Middleton's  (3a)  result,  the  power  spectral  density  of  the  random 
pulse  train  is 

oo  oo 

W(f)  =  w;;(n[  I  Ky(kT^)  e^'"^'^o  +  ^  I 

k=-CD  m=o 

(8) 

where 

W^(i)  =  average  power  spectral  density  of  the  individual  pulses 
Ky(t^-t2)  =  (yi“yi)(y2"y2^  ”  noise  autovariance  function 
y  =  average  value  of  the  pulse  amplitudes 
T^  =  sampling  period 
6  =  Dirac  delta  function. 

If  T^  is  large  compared  with  the  noise  fluctuation  time  -- 
roughly  the  mean  spacing  between  successive  zero-crossings  of 
the  noise  --,  then 

_  -  _2  °° 

W(f)  =  W^(f)  (y^  -  7^)  +  ^  'i  6(f  -  ^  )  . 

o  o 

m=o 


(9) 
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The  spectrum  consists  of  a  continuous  part  whose  distribution 
is  determined  solely  by  the  shapes  of  the  individual  pulses  and 
whose  intensity  is  proportional  to  the  variance  of  the  pulse 
amplitudes.  In  addition  there  is  a  discrete  spectrum  with  components 
at  multiples  of  the  sampling  rate  whose  intensity  is 

proportional  to  the  square  of  the  average  pulse  amplitude.  A 
typical  W(f)  is  shown  in  Figure  5.  If  y  =  0,  then  spikes  in  W(f) 
vanish.  For  the  white  noise  source  used  in  the  digital  system  the 
frequency  region  from  near  zero  to  a  little  less  than  i/T^^  cps 
was  extracted  by  a  passive  linear  low-pass  filter. 


2.  3 

0  0 

FREQUENCY 


Figure  5  The  Power  Spectral  Density  for  a 

Periodic  Pulse  Train  Having  Statis¬ 
tically  Independent  Pulse  Amplitudes, 
Pulse  Duration  t,  and  Period  T  . 
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Consider  the  random  pulse  train  x(t),  see  Figure  4,  as  the 
input  to  a  linear  filter  with  impulse  response  h(t)  as  is  shown  in 
Figure  6.  Let  x(t)  be  represented  by  a  weighted  sum  of  narrow  pulses; 
oo 

x(t)  =  i  ■  <to) 

n=  -  OO 


NARROW  PULSES  WEIGHTED  IMPULSE 

RESPONSES 

Figure  6  The  Linear  Filter. 


Then,  if  we  assume  that  the  duration  ofu^(t)is  short  compared  with 
the  duration  of  h(t),  u^(t)  i  6(t)  and  the  output  V{t)  of  the  filter  becomes 

00 

oo 

V(t)  =  J  x(T)h(t-T)dT  =  y^h(t-nT^)  .  (H) 

-oo 

n=-oo 

The  filter  output  V(t)  consists  of  the  weighted  sum  of  impulse 
responses  whose  epochs  are  at  integral  multiples  of  T^,  as  shown 
in  Figure  7. 

What  may  be  said  regarding  the  amplitude  distribution 
of  V(t)  ?  In  general  very  little,  for  this  is  related  to  the  more 
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Figure  7 


Randomly  Weighted  Filter  Impulse  Responses 
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general  problem  of  determining  the  output  amplitude  distribution 
of  a  linear  system  driven  by  a  non-Gaussian  random  process. 
However,  the  form  of  Equation  (11)  suggests  the  application  of  the 
central  limit  theorem  of  probability.  The  central  limit  theorem 
asserts  that  under  general  conditions  the  sum  of  a  large  number 
of  independent  random  variables  approaches  a  Gaussian  distribution. 
That  is  to  say,  if  n  independent  impulse  responses  produce 
non -negligible  effects  on  V(t^),  and  if  n  is  sufficiently  large,  then 
the  ensemble  of  values  V(t^)  will  be  Gaussianly  distributed.  Since 
the  distribution  is  tending  toward  a  Gaussian  distribution,  we  can 
perhaps  help  to  minimise  n  by  providing  Gaussianly  distributed 
pulse  amplitudes.  In  any  case  this  is  most  convenient,  for  we 
need  only  sample  the  amplitudes  of  a  Gaussian  noise  voltage 
whose  fluctuation  time  is  approximately  one -tenth  the  sampling 
period  in  order  to  sufficiently  approach  statistically  independent 
samples.  That  is,  if  the  sampling  rate  is  10  Kcps,  then  we  need 
only  sample  Gaussian  noise  having  an  approximately  uniform  power 
spectral  density  from  approximately  zero  frequency  to  100  Kcps. 

The  amplitude  probability  density  function,  p(x),  for  the 
random  pulse  train  x(t)  whose  amplitudes  are  Gaussianly  distributed 
with  zero  mean  is  given  by  ^ 


+  (1  ■  )6  (x) 

o 


e 


(12) 
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Also,  the  characteristic  function  f(v)  of  p(x)  is  given  by 

f(v)  =  J  e*^  dP(x)  =  m-  exp  [  '  "TT" 

-00  o  L  J  o 

(13) 

where  P(x)  is  the  probability  distribution  function  of  x.  Notice 
that  if  T  =  T^,  then  x(t)  is  a  boxcar  waveform  and  is  Gaussianly 
distributed.  Of  course,  the  output  of  the  linear  filter  under  these 
conditions  will  also  be  Gaussianly  distributed. 

In  order  to  apply  the  central  limit  theorem  to  the  situation 
under  discussion,  the  impulse  response  of  the  linear  filter 
must  have  a  significant  width,  or  memory,  extending  out  to  nT^. 

The  value  of  n  required  will  normally  be  of  the  order  of  ten.  For 
our  particular  case  the  value  of  n  required  to  regenerate  a 
Gaussian  distribution  was  experimentally  determined  to  be  of  the 
order  of  3  or  4.  Hence,  for  a  sampling  rate  of  10  Kcps  a  low- 
pass  filter  having  a  5  Kcps  cut-off  frequency  will  suffice  for  the 
linear  filter. 

Figure  8  shows  the  measured  amplitude  distributions  compared 
with  points  representing  a  Gaussian  distribution  for  the  case  of  a 
sampling  rate  of  10  Kcps  and  a  low-pass  filter  having  a  5  Kcps 
cut-off  frequency. 
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p(v) 


0  V 


Figure  8 


The  Measured  Amplitude  Distribution 
p(v)  Compared  with  Points  Represent¬ 
ing  a  Gaussian  Distribution. 
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C.  THE  POWER  SPECTRA  CONSIDERED 

The  results  presented  in  this  report  are  concerned  with 
the  following  power  spectral  densities; 


These  particular  power  spectral  densities  were  selected  for 
investigation  because  they  are  representative  of  the  spectra  that  occur 
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in  many  branches  of  science  and  technology,  and  because 
theoretical  approximations  for  the  corresponding  Pq(t)  have 
been  reported  (2,  8,  19). 

The  corresponding  system  functions  Hj^OZirf)  were  synthesized 
on  an  analog  computer  as  indicated  in  Figures  9  and  10.  The  system 
function  H^(j2irf)  corresponding  to  the  Gaussian  power  spectral 
density  W^(f)  was  synthesized  by  stagger  tuning  9  single  tuned 
filters.  In  this  manner  the  poles  of  H^(j2iTf)  were  located  at  the 
theoretical  (33)  pole  positions  corresponding  to  a  9  th  order 
approximation  to  the  unrealizable  Gaussian  power  spectral  density. 
Figure  il  shows  a  comparison  of  the  sinusoidal  gain  of  the  filter 
structure  as  a  function  of  frequency  with  the  theoretical  Gaussian 
function.  The  comparison  is  excellent. 

The  Wiener-Khinchin  theorem  gives  the  normalized 
autocorrelation  function  for  the  above  power  spectra  as: 

GO 

r  W  (f)  cos  Zitfrdf 

I  W^(f)  df 

°  n  =  0,  1,2,  3,4 


The  results  of  the  Integrations  yield: 


n=l 


(20) 


thesis  of  the  System  Functions  H  (jZirf)  for  n  = 


42 


Figure  10  The  Synthesis  of  the  System  Functions  ® 
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-U).  |T  1 

Pi(t)  =  (1  +  lu^  |T|)e 

(21) 

(22) 

r  1  1 

P2<t)  =  [1  +  UJ2ITI  +  — —  ]e 

P3(t)=  [i+U)j|t|  -  2(uj3t)^  +  3  (lUjI'i  |)^] 

,  -‘"3ITI 

e 

(2  3) 

P4(t)  =  cos  uu^T  e 

(24) 

The  quantity  ud^t  appearing  in  the  normalized  autocorrelation 

functions  is  defined  as  normalized  time  and  is  denoted  by  u  .  All 

n 

probabilities  and  probability  densities  reported  in  this  report  are 
plotted  with  respect  to  normalized  time  u^  = 

The  initiate  commands  discussed  in  Section  II-A  above  were 
spaced  approximately  one  second  apart.  Accordingly,  in  order  to 
obtain  samples  of  zero-crossing  intervals  which  are  sufficiently 
Independent,  the  time  constants  and  frequencies  associated  with 
the  power  spectral  densities  were  chosen  as; 


=  Tj  =  10m  sec 

(25) 

=  5  m  sec 

(26) 

=  80  cps 

(27) 

=  30.8  cps 

(28) 

Also,  for  experimental  convenience  the  synthesized  power 
spectral  densities  had  a  low  freqiency  cutoff  of  approximately  1  cps 
and  a  high  frequency  cutoff  of  8.5  kc. 
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D.  RANDOM  POINT- PROCESS  TRANSFORMATIONS 
Consider  the  linear  system  shown  in  Figure  12. 


Figure  12  The  Linear  System. 


The  power  spectral  density  for  the  Gaussian  random  process  x(t) 
is  given  by 

W^(f)  =  |H(j2xf)|^  .  (29) 

Similarly,  the  power  spectral  density  for  y(t)  is  given  by 
Wy(f)  =  |H(j2Trf)|^|Dj(j2Trf)|^  =  W^(f)  |  D^(j2nf)  |  ^ 

(30) 

If  D|(J2irf)  is  the  system  function  representing  an  ideal  differentiator 
then  the  stationary  points  of  x(t)  and  the  zero  points  of  y(t)  occur 
simultaneously . 

If  y(t)  is  a  Gaussian  process  having  arbitrary  power  spectral 
density  Wy(f),  the  average  number  of  zero  points  per  second  is 
given  by  (27) 
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(32) 


This  last  equation  agrees  with  Rice's  result  (27). 

The  average  number  of  inflection  points  per  second  for  a 
Gaussian  process  x(t)  having  power  spectral  density  ^^(^)  can  be 
determined  in  a  similar  manner.  Thus,  if  also  D2(j2iTf)  =  j2'iTf,  then 
the  average  number  of  inflection  points  per  second  for  a  Gaussian 
process  x(t)  having  power  spectral  density  W^(f)  must  be 


(33) 


The  same  transformation  can  bo  extended  to  higher  derivatives. 

A  method  for  synthesizing  a  Gaussian  process  y(t)  having 
zero  points  corresponding  in  time  with  the  stationary  points  of  a 


-47- 


Gaussian  proceas  x(t)  having  a  Batterworth  power  spectral 
density  is  shown  in  Figure  13. 


Figure  13  Transformation  of  Stationary  Points  to  Zero  Points. 


The  circuit  having  the  system  function  D(j2iTf)  is  shown  in 
Figure  14. 


5K  0.02»Jf 

• — WV— {(--■' 


lOOfipf 


IM 

f-WV-4 


o 


Figure  14  The  Synthesis  of  D(J2iTf). 
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The  Byetem  function  D(j2iTf)  is  given  by 

D(j2,rf)=  (34) 

[l+j  10“*u)] 

Accordingly,  D(J2iTf)  hae  the  property  of  a  differentiator  since 
|D(j2iTf)|  behaves  like  V"  to  well  within  ±  i  db  in  the  frequency 
range  of  interest  (0-500  cps). 

Taking  D(j2irf)  as  2jl0"^u)  we  find  that 


This  final  ratio  was  verified  satisfactorily  by  using  an  RMS 
voltmeter. 

Oscilloscope  photographs  illustrating  the  transformation  of 
stationary  points  to  zero  points  by  the  operator  D(j2iTf)  are  shown 
in  Figure  15.  The  probability  densities  associated  with  the  intervals 
defined  by  the  mathematical  stationary  points  of  a  random  process 
consisting  of  a  sine  wave  plus  a  Gaussian  process  were  determined 
experimentally  by  using  the  operator  0021^). 

Figure  16  shows  a  complete  block  diagram  of  the  measuring 


and  synthesis  systems. 


50 


Figure  16  Block  Diagram  of  the  Synthesis  and  Measuring  Systems. 
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IV.  MATHEMATICAL  DEDUCTIONS  FROM  THE 
PROBABILITY  DENSITIES 

A.  Kohlenberg  (4),  P.  I.  Kuznetsov,  P.  L.  Stratonovich  and 

V.  I.  Tikhonov  (5),  J.  A.  McFadden  (2),  and  perhaps  others,  have 
derived  some  basic  identities  relating  p(n,T),  the  probability  that 

a  given  interval  (t,  t  +  t)  contains  exactly  n  points,  and  PjjtT),  the 
probability  density  of  the  interval  between  the  mth  and  (m+n+l)th 
points.  These  identities  apply  to  general  stationary  point -proc e s se s . 
Using  McFadden's  general  result  (26)  we  have  for  a  point  set 
defined  by  the  zeros  of  a  random  process: 

oo 

p'(n.T)  =  -p  J  [  Pjf)  -  +  Pn.zW)] 

(37) 

where 

P  =  the  expected  number  of  zeros  per  unit  time 

P  (t)  =  the  probability  density  function  for  the  sum  of 
n+l  successive  zero-crossing  intervals 

p'(n.  t)  =  the  derivative  with  respect  to  t  of  the  probability 
that  a  given  interval  t  contains  exactly  n  zeros. 

Also,  by  definition,  we  have 

Z(n+1.T)=  J  P(l)di  (38) 

T 

where 

Z(n+1,T)  =  the  conditional  probability  that  the  (n  +  l)th  zero- 
crossing  point  from  a  given  zero-crossing  point 
in  dr  occurs  after  the  time  dr  +  t. 
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For  n  =  0,  and  the  initial  condition  p(0,  0)  =  i.  Equation  (37)  yields 

P(0,t)  =  P  dr^  J  ^  P^(l)  di  -T  +i  .  (39) 

o 

For  n  =  1  and  the  initial  condition  p(l,  0)  =  0,  Equation  (37)  yields 

'■  (40) 

For  n  =  0  Equation  (38)  becomes 

oo 

Z(1.t)=  f  Pjt)dl  .  (41) 

T 

For  n  =  1  Equation  (38)  becomes 
oo 

Z(2.T)=r^  Pl(i)dJ  (42) 

Equations  (39),  (40),  (41),  and  (42)  were  used  to  deduce 

p(0,  Uj^),  p(l,  u^),  Z(l,  u^)  and  Z(2,  u^)  for  n  =  1,2,3  and  4. 

Similar  expressions  hold  for  the  point  set  defined  by  the 

zero  points  of  a  random-process  consisting  of  a  sine  wave  signal 
f 

of  frequency  ^  plus  Gaussian  noise  such  that  the  signal-to-noise 
power  ratio  equals  a.  That  is: 

Pa(0*T)  =  Pa  /q '*■"1  J 

(43) 


I 
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s  . 


Pji.T)= 


dT.  f  P.(J.a)di-2  [’"dT,  P  (i.a)dJ+T 


00 

Z^(2,t)  =  J  P^(i.a)di 


(44) 

(45) 

(46) 


whera 


P  =  the  expected  number  of  zeros  per  unit  time  when 
^  the  signal  to  noise  power  ratio  equals  a. 

Equations  (43).  (44),  (45)  and  (46)  were  used  to  deduce 

p^(o,u^).  Pj^(l.u^).  Z^(l,u^)  and  Zj^(2,  u^)  for  a  =  0,  .2,  1  and  4. 

J 

Similar  expressions  also  hold  for  the  point  set  defined  by  the 

mathematical  stationary  points  of  a  random  process  consisting  of  a 

f 

sine  wave  signal  of  frequency  plus  Gaussian  noise  such  that  the 
signal -to-noise  power  ratio  equals  a.  That  is; 


m-(0,T)  =  y 


+  1 
(47) 


(.T  ,.T, 

.♦T  J 

Va 

dr,  'M,(f,a)dJ-2 

0  *  *  0  * 

dr,  J  (l,a)di+T 

•*  0  '  '*0  ° 

(48) 

SJl.T)  = 

&  «| 

f 

T 

M^(l.a)  di 

(49) 

s.(a,T)  =  J 

T 

Mj(f  ,a)  dj 

(50) 
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where 

a  the  expected  number  of  stationary  points  per  unit 
time  when  the  signal-to-noise  power  ratio  equals  a. 

Equations  (47),  (48),  (49).  and  (50)  were  used  to  deduce 
mj^(0,u^),  m^(l,u^),  S^(l,u^)  and  S^(2,u^)  for  a  s  0,  .2,  1  and  4. 
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V.  PROPERTIES  OF  THE  PROBABILITY  DENSITIES 
A.  EXPECTATIONS 

Consider  the  point  set  shown  in  Figure  17  and  defined  by 
the  zero-crossing  points  of  a  Gaussian  process.  Let  X,  Y  be  two 
successive  zero-crossing  intervals  and  let  Z  =  X  +  Y .  Prom  symmetry 
the  expectation  E(X)  equals  the  expectation  E(Y).  Hence,  E(Z)s:2E(X). 
Accordingly,  the  expectations,  associated  with  the  P^(u^) 

densities  are  related  to  the  expectations,  Eq(u^), associated  with 
the  Po^'^n^  densities  by: 

E^(u^)  =  2Eq(u^)  for  n  =  1. 2.  3,  4  (51) 

The  exact  theoretical  expectations,  Eq(u^),  associated  with 
the  densities,  or  the  average  values  of  successive  zero-crossing 

intervals,  are  given  by  (27) 


SuccsMtvs  lnt»rvals,t.  ,  Offintd  By  Adjoctnl  Stationary  Painti 


Figure  17  Point  Sets  Defined  by  the 
Zero-Crossing  Points  and 
the  Stationary  Points  of  a 
Gaussian  Process. 
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for  n  =  1, 2,  3,  4 


pVKu  )  is  the  second  derivative  o£  the  normalized  autocorrelation 
n  n 

function  The  results  for  n  =  1,2,3,  and  4  are: 


Eq(u^)  =  V  (53) 

Eq(u2)  =  5,4414  (54) 

Eq(u3)  =  1.4050  (55) 

Eq(u^)  =  3.  1341  (56) 


From  Equation  (51)  the  theoretical  expectations,  Ej(u  ), 

1  n 

associated  with  the  P^(u^)  densities  are  equal  to  2Eq(u^)  for 

n  =  1, 2,  3,  4.  The  expectations  Eq(u^)  and  E|(u^)  were  also  measured. 

They  are  compared  with  the  corresponding  theoretical  values  in  the 

tables  accompanying  Figures  18,  19,  22,  2  3,  26,  27,  30,  and  31.  The 

comparisons  are  satisfactory  in  all  cases. 

Similarly,  the  expectations,  E^(u^,  a), associated  with  the 

Pl(Uo'*^)  densities  are  related  to  the  expectations,  Ejj(u^,  a),  associated 

with  the  P  (u^,  a)  densities  by; 

o  o  ' 

Ei(Uo,a)  =  2Eo(Uo,  a)  for  all  a  .  (57) 

The  exact  theoretical  expectations,  Ejj(u^,  a),  associated  with 
the  a)  densities  follow  from  a  general  result  given  by  Rice  (34) 

and  Middleton  (3,  3a): 
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I-  oo 


oo 


for  all  a 


I  ^ 


(-baf 


Ln=o 


(n'T 


-1 

lF^(^;n+l;-a)j 


(58) 


where 


(q)„  =  (a)  (a  +  1)  (a  +  2)  .  .  .  (q  +  n  -  1) 
(a)o  =  1 

_  S  _  average  elne  wave  power 
*  “  TJ  “  average  noise  power 
3e 


Sin 


TT 


4Sin 


Tl 


=  the  confluent  hypergeometric  function. 
The  reeulte  for  a  =  0,  .2,  1  and  4  are: 


E  (u  ,  0)  3  5.2587 
o  o 

E  (u  ,  .2) =  5.4065 
o  o 

E^(u^,  1)  =  5.8286 

E  (u  ,4)  =  6.2506 
o  o 


(59) 

(60) 
(61) 
(62) 


For  small  values  of  "a",  the  random  process  consisting  of  a 

f 

sine  wave  signal  of  frequency  plus  Gaussian  noise  Is  approximately 

Gaussian.  Accordingly,  one  could  compute  easily  the  approximate 

theoretical  expectations,  a),  associated  with  the  Pq('1q>  a) 

densities  by  using  an  equation  analagous  to  Equation  (52).  That  is: 

-1/2  ,  ...  ,1/2 


E*(u^.  a)  =  X  [  -R"(0.  a)]  =  2it  /F  [  ] 


(63) 


R"(u  ,a)  is  the  second  derivative  of  the  normalized  autocorrelation 

°  fo 

function,  R(u  ,  a),  defined  by  the  power  spectral  density  W  (f)+S6(f— y). 

O  O  M 

R(u^,a)  is  given  by: 


-58- 


R(u^.a)  = 


u 

Po(u^)+aco8 -2^ 

1  +  a 


(64) 


Equation  (63)  agrees  with  Bendat's  result  (15). 
The  results  for  a  =  0,  .2.  1  and  4  are: 


E*(u  ,0)  =  5.2587 
o'  o 

(65) 

E*(u  ,  .2)  =  5.  3951 
o  o 

(66) 

E*(u^,  1)  =  5.7031 

(67) 

E*(u  ,  4)  =  6.  0  306 

0  o 

(68) 

From  Equation  (57)  the  exact  theoretical  expectations, 
E|(u^,a),  associated  with  the  P^(u^,a)  densities  are  equal  to 
2E^(Uo,  a)  for  all  a.  The  expectations  a)  and  Ej(Uj^,  a) 

for  a  =  1,  .2,  1,  and  4  were  also  measured.  They  are  compared 
with  the  corresponding  exact  theoretical  values  in  the  tables 
accompanying  Figures  34,  35  ,  38  ,  39  ,  42  ,  43,  46,  and  47.  The 
comparisons  are  satisfactory  in  all  cases. 

Finally,  the  expectations,  C^(u^,  a),  associated  with  the 
M^(u^,a)  densities  are  related  to  the  expectations,  ^QC^^^’a), 
associated  with  the  M^(u^,  a)  densities  by: 

6^(Uo'®^)  =  2eo(%»a)  all  a  (69) 

The  exact  theoretical  expectations,  ^^(^>^1  a),  associated  with 
the  M^(u^,  a)  densities  follow  by  applying  the  transformation 
associated  with  the  ideal  differentiator  operator  to  the  general 
result  given  by  Rice  and  Middleton: 
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oo 

£  (u  ,a)=  r  u  M  (u  ,  a)du 
o'  o  J  o  o'  o  '  o 
o 


Ln=0 


where 

c 


Sin 


5tt 

Ti 


1 6Sin 


for  all  a 


(70) 


The  results  for  a  =  0,  .2,  1  and  4  are: 


£  (u  ,  0) =  3.7765 
o  o 

(71) 

£^(u^,  .2)  =  3.9429 

(72) 

£  (u  ,  1)  =  4,5467 
o'  o 

(7  3) 

£  (u  ,  4)  =  5.8138 
o  o 

(74) 

Again,  for  small  values  of  "a"  the  random  process 

f 

consisting  of  a  sine  wave  signal  of  frequency  plus  Gaussian 
noise  is  approximately  Gaussian.  Accordingly,  one  could  compute 
easily  the  approximate  theoretical  expectations,  a),  associated 

with  the  M^(u^,a)  densities  by  using  an  equation  analagous  to 
Equation  (52).  That  is; 


eX,a)  =  ,[-R'J(0,,)] 


-1/2 


1  +  ab 
1  +  ac 


(75) 

R'J  (u^,  a)  is  the  second  derivative  of  the  normalized  autocorrelation 
function,  R  (u  ,a),  defined  by  the  power  spectral  density 

8  0  m 

2  2  ^ 

U)'  ‘'Wo(f)  +  U)  S6(f  ■  -j-  )•  This  latter  power  spectral  density  results  by 
applying  the  transformation  associated  with  the  ideal  differentiator 
operator  to  the  power  spectral  density 
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f 

W^(£)  +  S6(f  ^  ).  R^(u^,  a)  and  R(u^,  a)  are  related  by: 


R‘'(«  .a) 

R  (u  ,a)  = - for  all  a  .  (76) 

"  ^  R"{0.a) 

The  results  for  a  =  0,  .2,  1,  and  4  are: 

C*(V0)  =  3.776s  (77) 

e*(u^.  .2)=  3.9341  (78) 

6  *(vio.  1)  =  4.3994  (79) 

C*(Uo.  4)  =  5. 1910  (80) 


From  Equation  (69)  the  exact  theoretical  expectations, 

6  ^(u^,a),  associated  with  the  M|(vt^,  a)  densities  are  equal  to 
28  ^(u^,  a)  for  all  a.  The  expectations  8^(u^,  a)  and  8^(u^,a) 
for  a  =  0,  .2,  1,  and  4  were  also  measured.  They  are  compared 
with  the  corresponding  exact  theoretical  values  and  the  approximate 
theoretical  values  in  the  tables  accompanying  figures  50,  51,  54, 

55,  58,  59,  62,  and  63.  The  comparisons  are  satisfactory  in 
all  cases. 

If  one  represents  the  sine  wave  as  the  limit  of  a  narrow-band 

f 

Gaussian  process  centered  at  as  the  bandwidth  approaches  zero, 

then  it  becomes  clear  that  the  expectations  E  (u  ,  a)  and  8  (u  ,a) 

o  o  o  o 

are  exact  for  the  case  of  a  random  process  ^  (t,a)  of  the  form: 
e  (t,a)  =  R(t)  Sin  [  -^  +  0(t)  ]  +  n(t) 


(81) 
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In  this  last  equation  the  random  variable  has  a  Rayleigh 

probability  density  function,  6^  is  uniformly  distributed  in  the 

interval  (o,  Zit),  and  n(t)  is  Gaussian  noise  having  a  power  spectral 

density  W  (f).  In  other  words  £  (u  »a)  and  S  (u  ,  a)  are  exact  for 
'  o'  '  o  o  o  o 

the  Gaussian  process,  4  (t,  a),  defined  by  a  very  narrow-band 

f^ 

Gaussian  signal  centered  at  plus  Gaussian  noise  having  the 
power  spectral  density  W^(f).  Equation  (81)  should  be  compared 
with  random  process,  defined  by  a  sine  wave  signal 

plus  Gaussian  noise  which  has  the  form; 

4{t,a)  =  Sin  I  +  0(t)j  +  n(t)  where  is  a  constant. 

(82) 

B.  STANDARD  DEVIATIONS  AND  CORRELATION 
COEFFICIENTS 

Since  Z  =  X  +  Y,  the  variance  of  Z  is  given  by: 

Var(Z)  =  Var{X)  +  Var(Y)  +  2Cov(X,  Y)  (83) 

where 

Cov(X,Y)  =  E(XY)  -  E(X)  E(Y) 

The  correlation  coefficient,  K(X,  Y),  of  X,Y  is  given  by 


K(X,  Y)  = 


Cov(X,  Y) 
V^Var(X)Var(Y) 


From  symmetry,  we  have 
Var(X)  =  Var  (Y) 


(84) 


(85) 
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Hence, 

Var  (Z)  =  2  [  1  +  k]  Var  (X)  (86) 

Let  K  be  the  correlation  coefficient  between  two  successive 
n 

zero-crossing  intervals  of  a  Gaussian  process  having  the  power 

2 

spectral  density  W^(f).  Then,  the  variances,  o^(u^),  associated 
with  the  densities  are  related  to  the  variances,  o^(u^), 

associated  with  the  densities  by: 

0  j(Un)  =  2  1^  1  +  J  a^Cu^)  for  n  =  i,  2,  3,  4  (87) 

Similarly, 

o5(u^,0)  =  2  [l  +  K^]  oj(u^,0)  (88) 

d5(u^,  0)  =  2  [  1  +  ?Co]  D^(u^,  0)  (89) 

z  z  z 

In  the  latter  two  equations  a  j('»Q.  a),  ci  ^(u^.  a),  D^(u^,  a)  and 
D^(u^,  a)  are  the  variances  associated  with  the  respective  densities 
P^(Uq,  a),  P^(u^.  a),  M^(u^,  a)  and  M^(u^.  a).  Also  and  are 
the  correlation  coefficients  between  two  successive  zero-crossing 
intervals  of  the  Gaussian  processes  having  the  respective  power 
spectral  densities  W^(f)  and  lu^W^(f). 

McFadden  (32)  recently  derived  the  following  equations 
valid  for  a  general  random  process: 
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CD  O  j  ^ 

=  1  r(T)dT  =  ^  [  T  +  I  ]  t’O) 


EqCt) 


B=  J  [U(T)-  '  ]d,=  ^[*  +  £  P„]  -  J 


e;(t) 


where 

r(r)  =  the  normalized  autocorrelation  function  of  the 
infinitely  clipped  random  process 

=  variance  of  successive  zero-crossing  intervals 

■I 

Eo(t)  ~  ~  expectation  of  successive  zero-crossing  intervals 

p  a  correlation  coefficient  between  the  ith  and  (i+n)th 
”  zero-crossing  intervals 

U(T)dT  =  the  conditional  probability  that  a  zero  occurs 
between  t+r  and  t+r+dr  given  a  zero  at  time  t. 

If  one  assumes  that  successive  zero-crossing  intervals  form 

a  Markov  chain  in  the  wide  sense,  then  p  =  and 

n  1 


o 

TfTi 


for  all  I  p .  I  <1 


B  =  — ^  1  ).  (93) 

Equations  (92)  and  (93)  were  used  by  McFadden  (2)  to  compute 
and  p^  for  various  Gaussian  processes.  We  used  Equations  (92) 
and  (93)  to  compute  and  p^  for  most  of  the  Gaussian  processes 
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conaidcred  in  this  report.  In  some  cases  our  results  differ 
somewhat  from  those  reported  by  McFadden.  The  standard 
deviations  resulting  from  the  Markov  assumption  together  with 
the  standard  deviations  computed  from  Equations  (87),  (88)  and 
(89)  are  compared  with  the  experimental  standard  deviations  in  the 
tables  accompanying  Figures  18,  19,  22,  23,  26,  27,  34,  35,  50, 
and  51.  From  the  comparisons  of  the  standard  deviations  in  the 
tables  accompanying  Figures  26,  27,  50,  and  51,  we  conclude 
that  successive  zero-crossing  intervals  of  some  Gaussian  processes 
do  not  form  a  Markov  chain  in  the  wide  sense. 

Notice  that  the  Markov  assumption  Implies  that  if  >  0 
then  p,^  ^  ^  Also,  the  Markov  assumption  implies  the 

inequalities  given  by  Equations  92  and  93.  These  restrictions  are 
not  clearly  Implied  by  the  Gaussian  character  of  the  process  and 
depend  only  on  the  Markov  assumption.  Accordingly,  one  might 
surmise  that  the  successive  zero-crossing  intervals  of  some 
Gaussian  processes  do  not  form  a  Markov  chain  in  the  wide  sense. 

From  Equations  (90)  and  (91)  we  have  that; 


o 


(96) 
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By  applying  Equation  (96)  to  most  of  the  Gaussian  processes 
considered  in  this  report  and  by  using  the  experimental  variances 
for  successive  intervals,  the  experimental  correlation  coefficients 
between  two  successive  intervals  were  computed.  The  results 
are  given  in  the  tables  accompanying  Figures  18,  22,  26,  34,  and 
50. 

If  the  variances  for  the  sum  of  two  successive  intervals  are 
computed  from  Equations  87,  88,  and  89  by  using  the  experimental 
variances  for  successive  Intervals  and  the  experimental  correlation 
coefficients  between  two  successive  intervals, the  results  agree 
with  the  measured  variances  for  the  sum  of  two  successive 
intervals  to  within  the  experimental  error.  Hence  the  experimental 
values  of  the  correlation  coefficient*  between  two  successive 
intervals  yield  results  consistent  with  experiment.  This  serves 
as  partial  justification  for  the  assumption  given  by  Equation  (95). 

If  one  assumes  that  pj^  =  for  n  s:  2,  then  Equations 
(90),  (91),  and  (86),  can  be  used  to  compute  approximate 
theoretical  values  of  c^,  and  0|  for  Gaussian  processes.  For 
the  Gaussian  processes  considered  in  this  report,  these  approximate 
theoretical  values  compare  with  our  experimental  values  slightly  better 
than  the  corresponding  theoretical  values  resulting  from  the  Markov 
assumption. 
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For  a  Gaussian  process,  Rice  (i)  derived  U(T)dT,  see 

Appendix,  the  conditional  probability  that  a  Kero  occurs  between 

t  +  T  and  t  +  T  +  dr  given  a  zero  at  t.  Hence, 

00 

U(t)=  Y. 

n=0 

Also,  for  a  Gaussian  process.  Rice  (1)  derived  Q(T)dT,  see 
Appendix,  the  conditional  probability  that  a  downward  zero- 
crossing  occurs  between  t  +  t  and  t  +  t  +  dr  given  an  upward 
zero-crossing  at  t.  Hence 
00 

m=0 

and  ^ 

00 

U(t)-Q(t)=  I  P2„„(r)  (99) 

m=0 

From  Figure  3i  and  Equation  (98)  we  see  that  the  first 
portion  of  Q(u^)  generates  P^(u^).  Similarly,  from  Figure  32 
and  Equation  (99)  we  see  that  the  first  portion  of  U(u^)  -  Q(u^) 
generates  P^(u^).  The  first  normalized  portion  of  the  Q(u^) 
function  was  used  to  compute  approximate  theoretical  means  and 
variances  associated  with  the  P^(u^)  density.  Similarly, 
the  first  normalized  portion  of  the  U(u^)  -  Q(u^)  function  was  used 
to  compute  approximate  theoretical  means  and  variances  associated 
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with  the  density.  By  using  these  computed  variances  in 

Equation  (87),  the  approximate  theoretical  correlation  coefficient 
between  two  successive  intervals,  K^,  was  computed.  Because 
of  the  unfavorable  propagation  of  errors  associated  with 
Equation  (87),  should  be  regarded  as  a  rough  estimate. 

These  approximate  theoretical  results  are  compared  with 
other  results  in  the  table  accompanying  Figures  30  and  31. 

All  experimental  standard  deviations  were  computed  by 
using  the  theoretical  mean  values.  In  general,  these  standard 
deviations  would  decrease  somewhat  if  one  were  to  use  the 
experimental  mean  values  in  the  computations. 
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VI.  PROBABILITIES  AND  PROBABILITY  DENSITIES  DEFINED 
BY  THE  ZERO  POINTS  OF  GAUSSIAN  PROCESSES 
A.  POWER  SPECTRAL  DENSITY  W^(f) 

Probabilities  and  probability  densities  defined  by  the 
zero-crossing  points  of  a  Gaussian  process  having  a  power  spectral 
density  W|(f)  are  presented  in  Figures  18,  19,  20,  and  21.  Theore¬ 
tical  approximations  derived  by  S.  O.  Rice(l)  and  J.  A.  McFadden 
(2)  are  presented  in  Figures  18  and  19  for  comparison  purposes. 
Figure  18  also  presents  for  comparison  an  exponential  density 
which  was  suggested  by  R.  R.  Favreau,  H.  Low,  and  I.  Pfeffer  (19) 
as  the  exact  probability  density.  From  this  latter  comparison 
we  see  that  the  experimental  Pq(U|)  and  the  exponential  curve 
deviate  significantly  in  the  neighborhoods  of  u^  s  1  and  u^  =  5. 
Accordingly,  we  conclude  that  the  exact  probability  density  is  not 
the  suggested  exponential. 

From  Equation  (98)  we  have  that; 

Q(t)  i  Pq(t)  for  all  t  (100) 

Accordingly,  Rice's  Q(U|)  function  serves  as  an  upper  bound 

for  Pq(U|)  And  is  presented  in  Figure  18.  Recently,  M.  S.  Longuet- 

Higgins  (8)  deduced  bounds  for  Pq(0)  theoretically  as  0.  382  < 

Pq(0)  <  Q(0)  =  0.  406  and  deduced  the  approximate  value  Pq(0)  =  0.  385. 
The  inequality  0.  382  <  Pq(0)  was  used  by  Longuet- Higgins  (8)  to 
disprove  rigorously  the  exponential  conjecture  of  Favreau  et  al. 

From  Equation  (99)  we  have  that: 


Figure  20  p{n,  u^)  is  the  Probability  that  a  Given  Interval  u.  Contains  Exactly  n  Zero- 
Crossing  Points,  pfn,  U|)  and  White's  W_^n,  u.)  Probabilities  are  Compared 
for  the  Case  of  Gaussian  Noise  Having  a  Power  Spectral  Density 
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U(t)  -  Q(t)  i  P^(t)  for  all  t.  (101) 

Accordingly,  Rice's  U(u^)  -  Q(u^)  function  serves  as  an 
upper  bound  for  and  is  presented  in  Figure  19.  Longuet- 

Higgins  (8)  deduced  bounds  for  P^(0)  theoretically  as  0.  0645 
<  Pj(0)  <  U(0)  -  Q(0)  T  0.0726  and  deduced  the  approximate  value 
Pj(0)  i  0.0653. 

G.  M.  White  (18)  experimentally  detennined  the  following 
probabilities  concerning  the  zero-crossings  of  a  Gaussian  process 
having  a  power  spectral  density  W^(f): 

(1)  W^(0,Uj),  the  probability  of  having  no  zero- 
crossings  with  positive  slope  in  a  time  interval  u^. 

(2)  W^(l,Uj),  the  probability  of  having  one  zero-crossing 
with  positive  slope  in  a  time  interval  u^. 

In  general 

W^(0,t)  =  p(0,t) +y  p(1.t)  (102) 

and 

W^(1,t)  =  |•p(l,T)  +  p(2,t)  +  ^p(3,  t)  (103) 

or 

p(2,t)  s  W^(1,t)  -  (104) 

Equation  (102)  states  that  the  interval  can  contain  no  upcrossings 
in  two  ways:  Either  it  contains  no  crossings  at  all,  or  it  contains 
one.  If  it  contains  one,  then  with  conditional  probability  j  it  is 
either  upward  or  downward. 


-74- 


Equation  (103)  atatea  that  the  interval  can  contain  exactly 
one  upcroaalng  in  three  waya:  Either  it  containa  one  croaaing,  or 
two,  or  three.  If  it  containa  one,  then  with  conditional  probability 
it  ia  upward  or  downward.  If  it  containa  two,  then  one  of  them 
muat  be  upward.  If  it  containa  three,  then  one  or  two  of  them 
will  be  upward,  each  with  probability  j  . 

Figure  20  comparea  p(0,  u^)  -f  ^  p(l,  u^)  with  W^(0,  u^).  The 

compariaon  ia  excellent.  Alao,  Figure  20  preaenta  both  W^(l,  u^) 

and  p(l,U|).  Accordingly,  an  upper  bound  for  p(2,U|)  can  be 

readily  computed  from  inequality  (104).  Thla  upper  bound  agreea 

aatiafactorily  with  the  approximate  theoretical  initial  behavior 
0  818  2 

p(2,u^)~  — —  (uj)  Mdiich  waa  deduced  by  Longuet-Hlggina (8). 

Figure  19  comparea  the  convolution  of  Pq(uj)  with  itaelf, 
Po^'^1^  *  Thia  compariaon  aervea  to 

demonatrate  that  two  aucceaaive  aero-croaaing  intervala  are 
atatiatically  dependent.  A  theorem  ia  proved  in  Section  IX  which 
aaaerta  that  for  Gauaaian  proceaaea  the  aum  of  m  -f  1  aucceaaive 
zero-croaaing  Intervala  and  the  aum  of  the  next  n  -f  1  aucceaaive 
zero-croaaing  intervala  are  atatiatically  dependent.  The 
theorem  appllea  for  all  nonnegative  integral  m  and  n. 

B.  POWER  SPECTRAL  DENSITY  W2(f) 

Probabllitiea  and  probability  denaitiea  defined  by  the 
zero-croaaing  pointa  of  a  Gauaaian  proceaa  having  a  power  spectral 
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density  W2U)  are  presented  in  Figures  22,  23,  24,  and  25. 

Theoretical  approximations  derived  by  Rice  (1)  and  McFadden  (2) 
are  presented  in  Figure  22  for  comparison  purposes.  Figure  22 
also  presents  an  experimental  probability  density  reported  by 
Favreau,  et  al  (19).  Rice's  U(u2)-Q(u2)  function  is  presented 
in  Figure  23  for  comparison  purposes.  Figure  23  also  compares 
the  convolution  of  with  itself,  Po^“2^  * 

Pl(u2).  This  latter  comparison  was  intended  to  exhibit  the 
statistical  dependence  of  two  successive  intervals.  In  this  case 
the  convolution  function  is  nearly  equal  to  P^(u2).  However,  as  the 
theorem  of  Section  IX  asserts,  two  successive  zero-crossing 
intervals  are  still  statistically  dependent. 

The  theoretical  initial  values  of  Pq(u2)  and  P^(u2)  are 
zero  since  Rice's  U(u2)  function  is  initially  zero. 

C.  POWER  SPECTRAL  DENSITY  Wj(f) 

Probabilities  and  probability  densities  defined  by  the  zero¬ 
crossing  points  of  a  Gaussian  process  having  a  power  spectral  density 
W^(f)  are  presented  in  Figures  26,  27,  28,  and  29.  Theoretical 
approximations  derived  by  Rice  (1)  and  McFadden  (2)  are 
presented  in  Figure  26  for  comparison  purposes.  Figure  26  also 
presents  an  experimental  probability  density  reported  by 
Favreau,  et  al.  (19). 


Probability  Density  Function  for  the  Sum  of  Two  Successive  Zero 
Crossing  Intervals  of  Gaussian  Noise  Having  a  Power  Spectral  Density  W2(f). 

Compared  with  the  EIxperimental  Result  to  Exhibit  the 
Practical  Independence  of  Two  Successive  Intervals. 
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Figure  24  p(n,  U2^)  is  the  Probability  that  a  Given  Interval  U2  Contains  Exactly  n  Zero- 
Crossing  Points.  p(0,  u^)  and  p(l,  u  )  are  Compared  for  the  Case  of 
Gaussian  Noise  Having  a  Power  Spenral  Density 


Figure  25  Z(n,  u^)  i»  the  Conditional  Probability  that  the  nth  Zero-Crossing  Point  from 
a  Given  Zero-Crossing  Point  in  du.  Occurs  After  the  Time  du^+u^.  ZCl.u^) 
and  Z(2,  xi^)  are  Compared  for  the  Case  of  Gaussian  Noise  Having  a  Power 
Spectral  I^nsity 


81 


-82- 


Figure  28  p(n,  u^)  is  the  Probability  that  a  Given  interval  u-  Contains  Fxactly  n  Zero 
Cros^ug  Points.  p(0,  u,)  and  p(l,u,)  are  Compaired  for  the  Case  of 
Gaussian  Noise  Having  a  Power  Spenral  Density  W^(f). 


00 

d 

d 

Figure  29  Z(n,  u^)  is  the  ConditioDai  Probability’  that  the  nth  Zero-Crossing  Point 
from  a  Given  Zero-Crossing  Point  in  du~  Occurs  After  the  Time  du_+u 
Z(l,  U-)  and  Z(2,  u.)  are  Compared  for  tlie  Case  of  Gaussian  Noise 
Ha'ving  a  Power  S^Sectral  Density  W-(f). 
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Rlce'a  Q(u^)  function  serves  as  an  upper  bound  for  Pq(u^) 
and  is  presented  in  Figure  26.  Longuet-Higgins  (8)  deduced 
bounds  for  Pq{0)  theoretically  as  0.  612  <  Pq(0)  <  Q(0)  =  0.650 
and  deduced  the  approximate  value  Fq(0)  =  0.616.  Also,  Longuet- 
Higgins  (8)  deduced  bounds  for  P|(0)  theoretically  as  0.  103  < 

P^(0)  <  U(0)  -  Q(0)  =  0.  116  and  deduced  the  approximate  value 
P^(0) i  0.104. 

The  experimental  probabilities  p(0,  u^)  for  n  =  1, 2,  3 
satisfy  theorems  1,  3  and  Equation  (22)  of  D.  Slepian’s  (13)  recent 
paper  regarding  bounds  for  these  probabilities. 

Rice's  ^U(u^)-Q(uj)  j  function  is  presented  in  Figure  27  for 
comparison  purposes.  Figure  27  also  compares  the  convolution  of 
Pq(Uj)  with  Itself,  Pq(u2)  *  Pq(u^),  and  P^(uj).  This  latter 
comparison  serves  to  demonstrate  that  two  successive  zero- 
crossing  Intervals  are  statistically  dependent. 

The  experimental  probability  densities  Pq(u^)  for  n  =  1,2,3 
agree  very  well  with  theoretical  approximations  recently  published 
by  Longuet-Higgins  (8).  The  approximate  theoretical  initial 
values  of  Longuet-Higgins  were  used  in  the  normalization  of  the 
experimental  densities  ~  1> 

D.  POWER  SPECTRAL  DENSITY  W^(f) 

Probabilities  and  probability  densities  defined  by  the 
zero-crossing  points  of  a  Gaussian  process  having  a  power  spectral 
density  W^(f)  are  presented  in  Figures  30,  31,  32,  and  33.  Theoretical 
approximations  derived  by  Rice  (1 )  are  presented  in  Figures  30  and  31  for 
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U4  «  2'r 

Figure  31  P^(u^)  is  the  Probability  Density  Function  for  the  Sum  of  Two  Successive 

Zero-Crossing  Intervals  of  Gaussian  Noise  Having  a  Power  Spectral 
Density  WAf).  Q(u^)*Q(u^)  is  Compared  with  U(u^)-Q(u^)  to  Esdiibit 
the  Statistical  Dependence  of  Two  Successive  Integrals. 
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compariBonpurpoaea.  Figure  31  alao  comparea  the  convolution  of  the 
firat  normalized  portion  of  Q(u^)  with  itaelf,  Q(u^)*Q(u^).  and  the 
firat  normalized  portion  of  U(u^)  -  Q(u^),  In  accordance  with 
Equationa  (98)  and  (99),  thia  compariaon  aervea  to  demonstrate 
that  two  Bucceaaive  zero-croaaing  intervala  are  statistically 
dependent.  The  theoretical  initial  values  of  Pq(“4)  and  P|(u^) 
are  zero  since  Rice's  0(04)  function  is  initially  zero. 

The  results  in  Figures  32  and  33  were  computed  by  using 
Equationa  (39),  (40),  (41),  and  (42).  The  first  normalized  portions 
of  0(04)  and  11(04)  -  Q(u^)  served  as  Pq(u4)  and  Pi(n^) 
respectively. 

In  order  to  obtain  good  experimental  estimates  of  ^^(u^)  and 
0^(04),  we  see  from  Figures  30  and  31  that  the  experimental  densities 
P^(u^)  and  P|(n4)  must  be  measured  with  finer  resolution. 
Accordingly,  our  present  experimental  estimates  of  0  ^(u^)  and 
Oj(u^)  are  not  worth  reporting. 

Notice  that  a  "Q  factor*  of  6. 1  as  defined  in  Equation  (18)  is 
approximately  the  minimum  "Q  factor"  representing  a  "narrow -band 
system*  as  defined  in  communication  theory.  Accordingly,  the 
problem  of  theoretically  determining  Pq(t)  and  Pj(t)  for  narrow- 
band  Gaussian  processes  is  practically  solved  by  using  Rice's 
Q(t)  and^U(T)  -  Q(T)j functions  as  the  generating  functions  of  Pq(t) 
and  P|(t)  In  accordance  with  Equations  (98)  and  (99). 
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A  portion  of  the  above  work  dealing  with  Gaussian  processes 
has  been  published  recently  in  the  IRE  Transactions  on 
Information  Theory  (35). 
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Vn.  PROBABIUTIES  AND  PROBABIUTY  DENSITIES  DEFINED 
BY  THE  ZERO  POINTS  OF  A  SINE  WAVE  PLUS  A 
GAUSSIAN  PROCESS 

Probabilities  and  probability  densities  defined  by  the 

zero-crossing  points  of  the  random  process,  |(u,  a),  consisting 

f^ 

of  a  sine  wave  of  frequency  plus  Gaussian  noise  having  a 
power  spectral  density  W^(f)  are  presented  in  Figures  34  thru 
49.  The  ratio  of  the  average  sine  wave  power,  S,  to  the  average 
noise  power,  N,  is  denoted  by  the  parameter  a.  For  small  values 

g 

of  a  s  ,  the  random  process  |(u,  a)  is  approximately  Gaussian 
with  normalized  autocorrelation  function,  R(u^>  a),  given  by 
Equation  (64).  Accordingly,  Rice's  probability  function  Q  (u^,a)du^ 
approximates  the  conditional  probability  that  a  downward  zero¬ 
crossing  of  4(u,  a)  occurs  between  u  +  u^  and  u  +  u  +  given  an 

upward  zero-crossing  of  4(u,  a)  at  time  u.  In  accordance  with 
Equation  (98),  Q  (u^,  a)  serve  as  approximations  to  the  initial 
behavior  of  P^(u^,  a)  for  small  "a"  and  are  presented  in  Figures 
34,  38  ,  42,  and  46.  Figure  34  also  presents  McFadden's 
theoretical  approximation  for  Pq('1q»  0).  Q  (^Q.a)  du^  is  the 
exact  conditional  probability  that  a  downward  zero-crossing  of 
^  (u,  a)  occurs  between  u  -f  Uq  and  u  +  Uq  -i-  duQ,  given  an  upward 

l|l  ^ 

zero-crossing  of  ^  (u,  a)  at  time  u.  The  process  ^  (u,a)  is 
defined  by  Equation  (81)  and  the  process  ^(u,  a)  is  defined  by 
Equation  (82).  When  a  =  0  the  processes  |(u,  0)  and  ^  (u,  0)  are 


identical. 


Figure  39  P^(u^,  a)  is  the  Probability  Density  Function  for  the  Sum  of  Two  Successive  Zero- 

Crossing  Intervals  of  a  Random  Process  Consisting  of  a  Sine  Wave  of  Frequency- 
Plus  Gaussian  Noise  Having  the  Power  Spectral  Density-  W  (f).  The  Ratio  to  the 
Average  Sine  Wave  Power  to  the  Average  Noise  Power  is  fienoted  by  the  Parameter 
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Figure  44  “ft)  Probability  that  a  Given 

Imerv^  u  Contains  Exactly  n  Zero- 
Crossing  Points  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency 
f  /Z  Plus  Gaussian  Noise  Having  the 
Power  Spectral  Density  W  (f).  The  Ratio 
of  the  Average  Sine  Wave  Power  to  the 
Average  Noise  Power  is  Denoted  by  the 
Parameter  a. 
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Figure  48  Probability  that  a  Given 

Interv^  u  Contain*  Exactly  n  Zero- 
Crossing  Points  of  a  Random  Process 
Consisting  of  a  Sine  Wave  of  Frequency 
f  /Z  Plus  Gaussian  Noise  Having  the 
Power  Spectral  Density  W  (f).  The 
Ratio  of  the  Average  Sine  Wave  Power 
to  the  Average  Noise  Power  is  Denoted 
by  the  Parameter  a. 


c 

Similarly,  for  small  values  of  a  =  Rice's  probability 
41 

function  U  (uQ,a)duQ  approximates  the  conditional  probability  that  a 
zero  of  ^(u,  a)  occurs  between  u+Uj^  and  u+Ug+du^  given  a  zero  of 
^(u,  a)  at  time  u.  In  accordance  with  Equation  (99),  U  (u^,  a)-Q  (uQ,a) 
serve  as  approximations  to  the  initial  behavior  of  P^(u^,  a)  for  small  "a" 
andarepresentedinFigures  35,  39,  43,  and  47.  Notice  that  U  (uQ,a)du^ 
is  the  exact  conditional  probability  that  a  zero  of  ^  (u,  a)  occurs  between 
u-I-Uq  and  u+UQ+du^  given  a  zero  of  |  (u,  a)  at  time  u. 

Figure  35,  39,  43,  and  47  also  compare  the  convolution 
of  PQ(uQ,a)  with  itself,  Pq(uq,  a)*PQ(uQ,  a),  and  P^(uQ,a)  for 
a  s  0,  .2,  1,  4.  These  latter  comparisons  serve  to  demonstrate 
that  successive  zero-crossing  intervals  of  ^(u,  a)  are  statistically 
dependent . 

For  a  general  random  process,  McFadden  (28)  showed  that: 

'  I  (-<)“P„(T)  (105) 

n=0 

r"(T)  is  the  second  derivative  of  the  normalized  autocorrelation 

function  of  the  infinitely  clipped  random  process.  The  normalized 

autocorrelation  function  ^(uq,  a),  see  Appendix,  for  the  infinitely 

clipped  process  T  (u,  a)  can  be  computed  from  a  general  result 

given  by  Davenport  (40)  and  Middleton  (3a). 

Figure  46  compares  the  experimental  Pq('^qi  4)  with  a  first 
^11/^  4) 

portion  of - — -  Eq(Uq,  4)  in  accordance  with  Equation  (105). 
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The  comparison  is  excellent.  Figure  47  compares  the  experimental 

r"(uQ,  4) 

P^(uq,  4)  with  a  second  portion  of - ^ -  Eq(uq,  4)  in 

accordance  with  Equation  (105).  Again,  the  comparison  is 
excellent. 
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Vin.  PROBABILITIES  AND  PROBABILITY  DENSITIES  DEFINED 
BY  THE  STATIONARY  POINTS  OF  A  SINE  WAVE  PLUS 
A  GAUSSIAN  PROCESS 

Probabilities  and  probability  denaities  defined  by  the 
mathematical  stationary  points  of  the  random  process  ^(u,  a), 
defined  by  Equation  (82),  are  presented  in  Figures  50  thru  65.  For 

g 

small  values  of  a  =  the  rsmdom  process 
d^(u.a) 

tiK  a)= — -  (i06) 

is  approximately  Gaussian  with  normalized  autocorrelation  function, 
R^(uQ,a),  given  by  Equation  (76).  Accordingly,  Rice's  probability 
function  Q*(uQ,a)du^  approximates  the  conditional  probability  that 
a  downward  zero-crossing  of  ti(u,  a)  occurs  between  u  +  u^  and 
u  +  Ug  +  dug,  given  an  upward  zero-crossing  of  ti(u,  a)  at  time  u. 

In  accordance  with  Equation  (98),  Q^(ug,a)  serve  as  approximations 
to  the  initial  behavior  of  MQ(Ug,  a)  for  small  "a”  and  are  presented 

i|| 

in  Figures  50,  54,  58,  62.  Notice  that  Qj(ug.  a)  dUg  is  the  exact 
conditional  probability  that  a  downward  zero-crossing  of 

(i07) 

occurs  between  vi  +  Ug  and  u  +  Ug  +  dUg,  given  an  upward  zero- 

s 

crossing  of  q  (u,  a)  at  time  u. 

Similarly,  for  small  values  of  a  =  Rice's  probability 
function  U^(ug,  a)  dUg  approximates  the  initial  behavior  of  tho 


I 
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Figure  53  m  (n,  u  )  la  the  Probability  that  a  Given 
InfervaP  Contains  Exactly  n  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  f_/2  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  W  (f).  The  Ratio  of  the  Average 
Sine  Wave°Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 
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Figure  54  M  (u^,  a)  is  the  Probability  Density  Function  for  Successive  Intervals  Defined  by 

Adjacent  Stationary  Points  of  a  Random  Process  Consisting  of  a  Sine  Wave  of  Freq¬ 
uency  f  /2  Plus  Gaussian  Noise  Having  the  Power  Spectral  Densitv  W  (f).  The  Ratio 
of  The  Average  Sine  Wave  Power  to  the  Average  Noise  Power  is  r)eno*?ed  by  the 
Parameter  a.  Theoretical  Approximations  are  Compared  with  Ebcperimental  Results. 
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Figure  56  m  (n,  u  )  is  the  Probability  That  a  Given 
Infervaf  u  Contains  Exactly  n  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  fg/Z  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  W  W).  The  Ratio  of  the  Average 
Sine  Wave  r*ower  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 
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Figure  59M^(u^,  a)  is  the  Probability  Density  Function  for  the  Sum  of  Two  Successive  Intervals 

Defined  by  Adiacent  Stationary  Points  of  a  Random  Process  Consisting  of  a  Sine  Wave  of 
Frequency  f^/z  Plus  Gaussian  Noise  Having  the  Power  Spectral  Density  W  (f).  The  Ratio 
of  the  Average  Sine  Wave  Power  to  the  Average  Noise  Power  is  Denoted  by  the  Parameter  a. 
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Flgure  60  m  (n,  u  )  la  the  Probability  that  a  Given 
Interval  Containa  Ebcactly  n  Stationary 
Pointa  of  a  Random  Proceaa  Conaiating  of 
a  Sine  Wave  of  Frequency  f  /Z  Plua  Gauaa- 
ian  Noiae  Having  the  Power°Spectral  Denaity 
W^(f).  The  Ratio  of  the  Average  Sine  Wave 
Power  to  the  Average  Noiae  Power  ia  Denoted 
by  the  Parameter  a. 
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Figure  61  A  Random  Process  Consisting  of  a  Sine  Wave  of  Frequency  f  ,  2  Plus  Gaussian  Noise 
living  the  Power  Spectral  Density  W  (f)  Defines  S  (n,  u  ).  tl?e  Conditional  Probabili: 
tlmt  the  nth  Stationary  Point  from  a  CTven  Stationary  PoSit  in  du  Occurs  After  the 
Time  du^+u^.  The  Ratio  of  the  Average  Sine  Wave  Power  to  the'^A^'erage  Noise 
Power  IS  Denoted  by  the  Parameter  a. 
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Figure  64  u  )  i*  the  Probability  that  a  Given 

Inferval  Contains  Exactly  n  Stationary 
Points  of  a  Random  Process  Consisting 
of  a  Sine  Wave  of  Frequency  f  /Z  Plus 
Gaussian  Noise  Having  the  Power  Spectral 
Density  W  (f).  The  Ratio  of  the  Average 
Sine  Wave  Power  to  the  Average  Noise 
Power  is  Denoted  by  the  Parameter  a. 
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conditional  probability  that  a  stationary  point  of  ^(u,a)  occurs 

between  u  +  u^  and  u  +  Uq  +  du^  given  a  stationary  point  of  ^(u,  a) 

at  time  u.  In  accordance  with  Equation  (99),  U^iu^,  a)-Q^(uQ,  a) 

serve  as  approximations  to  the  initial  behavior  of  M^(uq,  a)  for 

small  "a"  and  are  presented  in  Figures  51,  55,  59,  and  63.  Notice 

that  U^(uQ,a)du^  is  the  exact  conditional  probability  that  a  stationary 
4t 

point  of  4  a)  occurs  between  u  )  and  u  +  Uq  +  dtiQ  given  a 

stationary  point  of  ^  (u,  a)  at  time  u. 

Figures  51,  55,  59,  and  6  3  also  compare  the  convolution 

of  M^(u^,  a)  with  itself,  M^(uQ,a)  *  M^(uq,  a),  and  M|(uQ,a) 

for  a  =  0,  .2,  1,  4.  These  comparisons  serve  to  demonstrate 

that  two  successive  intervals  defined  by  the  adjacent  stationary 

points  of  a)  are  statistically  dependent. 

The  normalized  autocorrelation  function  r^(Up,  a),  see 

Appendix,  for  the  infinitely  clipped  process  ti(u,  a)  can  be  computed 

by  applying  the  transformation  associated  with  the  ideal  differentiator 

operator  to  the  general  result  given  by  Davenport  and  Middleton. 

Figure  62  compares  the  experimental  M  (u,.,  4)  with  a  first 
r'^(Uo,a) 

portion  of - ^ 4)  in  accordance  with  Equation  (105). 

The  comparison  is  satisfactory. 

Figure  63  compares  the  experimental  M,(u„,  4)  with  a  second 

r"(ujj.a)  * 

portion  of - ^ accordance  with  Equation  (105). 

Again,  the  comparison  is  satisfactory. 


-128- 


IX.  STATISTICAL  DEPENDENCE  OF  ZERO- CROSSING 
INTERVALS 

The  recent  work  of  Longuet -Higgins  (8)  leads  one  to  the 
following  important  theorem. 

Theorem:  If  ^(t)  is  a  Gaussian  process  having  a  finite 
expected  number  of  zeros  per  unit  time,  then  the  sum  of  (m+i) 
successive  zero-crossing  intervals  of  ^(t)  and  the  sum  of  the 
next  (n+1)  successive  zero-crossing  Intervals  |(t)  are  statistically 
dependent  for  all  nonnegative  integral  m,n. 

Proof : 

For  Gaussian  processes  having  power  spectral  densities 

_4 

with  asymptotic  behavior  f  (the  singular  case),  Longuet -Higgins 
showed  that  F|(T)t  the  probability  density  function  for  the  sum  of 
(i  +  1)  successive  zero-crossing  intervals,  tends  to  a  positive 
value  as  T-*0.  Accordingly, 

)dl  =  0 

T-»0  T  -*0  T-^0  0 

(108) 

This  completes  the  proof  of  the  theorem  for  the  singular  case. 

For  Gaussian  processes  having  power  spectral  densities 

_4 

with  asymptotic  behavior  other  than  f  (the  regular  case),  Longuet- 

Higgins  showed  that 

1  (m+2)(m+3)-2 

Pm(T)  =  O  (t  )  as  t  -  0 


(109) 
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Hence,  as  t  0 

■i  (m+n+3)(m+n+4)-2 

and  .  . 

(n42){n4  3)+  —  (j-n+2)(m+3)-  3 
P„(T)  •  P„(T)  =  O  (/  ^  1  . 

(IH) 

Suppose  Equations  (110)and{lH)  were  equal  for  some  nonnegative 
integral  m,n.  Then. 

w(m+n+3)(m+n+4)-2  =  i(n+2)(n+3)f  ■i(m+2)(m+ 3)- 3 

(112) 

and, 

mn  +  m  +  n+  l  =  0  .  (113) 

Equation  (113)  is  clearly  false  for  all  noxuiegative  integral  m,n. 
Accordingly,  as  t  *♦  0 

This  completes  the  proof  of  the  theorem.  For  m  =  n  -  0,  the 
content  of  the  theorem  was  given  by  Palmer  (11)  and  McFadden 
(28). 

For  certain  Gaussian  processes,  the  convolutions  shown  in 
Figures  19,  27,  31,  35,  and  51  verify  the  truth  of  the  theorem 
for  m  =  n  =  0. 

In  order  to  observe  the  statistical  dependence  between  the  ith 
and  (l'fn)th  zero-crossing  intervals  of  a  random  process,  two  of 
the  experimental  systems  described  in  Section  HA  were  used  to 
display  5000  random  samples  of  the  ith  and  (i+n)th  zero-crossing 
intervals  on  an  oscilloscope.  The  skipping  discussed  in  Section 
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1 

IIA  was  set  at  11  or  12  each  with  probability  .  The  resulting 

intensity  pattern  on  the  oscilloscope  represents  an  approximation  to  ^ 

the  joint  probability  density  of  these  intervals.  If  the  Intensity 

pattern  corresponding  to  n  -♦oo  is  similar  to  the  intensity  pattern 

corresponding  to  n  =  n^,  then  one  can  conclude  that  the  ith  and 

(l+n^)th  zero-crossing  intervals  of  a  finite  memory  process  are 

practically  independent. 

For  a  Gaussian  process  havdng  power  spectral  density  W^(f), 

Figure  66  illustrates  the  statistical  dependence  between  the  ith 
and  (i+n)th  zero-crossing  intervals  of  the  process  for  n  =  1,2,  3,  4, 
and  00.  n  =  oo  represents  a  condition  such  that  the  two  zero- 

crossing  intervals  are  separated  by  approximately  71  average  « 

zero-crossing  intervals.  This  condition  was  convenient  experimentally 

and  corresponds  to  a  one  second  time  separation  between  the 

zero-crossing  intervals.  Notice  the  symmetry  of  the  intensity 

patterns  about  the  straight  line  passing  through  the  origin  with 

a  slope  of  45°.  This  symmetry  results  from  the  fact  that  the 

joint  probability  density  of  the  intervals  is  invariant  under  a  time 

reversal  of  the  random  process. 

The  conditional  mean  of  the  (i+n)th  zero-crossing  interval 
given  the  ith  zero-crossing  interval  is  found  from  Figure  66  by 
considering  narrow  vertical  strips  of  the  intensity  patterns  and 
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Figure  66:  Each  Intensity  Pattern  Represents  an 

Approximation  to  the  Joint  Probability 
Density  of  the  ith  and  (i+n)th  Zero-  ^ 
Crossing  Intervals  of  a  Gaussian  Process 
Having  Power  Spectral  Density  W^(f). 
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determining  the  center  of  gravity  for  each  strip  when  unit 
intensity  is  assigned  unit  mass.  A  line  joining  these  points 
is  called  a  regression  curve.  By  using  the  principle  of  least 
squares  one  can  determine  the  best  fitting  straight  line  to  the 
regression  curve.  The  slope  of  this  best  fitting  straight  line 
represents  the  correlation  coefficient  p^.  From  Figure  66 
we  see  that  when  n  2  3  the  regression  curve  is  approximately  a 
line  parallel  to  the  axis  of  abscissas,  and  the  conditional  mean 
value  is  approximately  equal  to  the  mean  value  of  the  zero- 
crossing  intervals.  Thus  when  n  a:  3,  =  0  or  the  zero- 

crossing  intervals  are  approximately  uncorrelated. 

By  comparing  the  intensity  pattern  for  n  =  00  to  the  others 
we  see  that  when  n  ^  4  the  statistical  dependence  practically 
vanishes.  However,  when  we  compare  the  intensity  pattern 
for  n  =  00  to  that  corresponding  to  n  s  1  we  see  significant 
dependence.  This  dependence  is  in  agreement  with  the  theorem 
of  this  section.  These  observations  may  be  useful  to  guide  the 
development  of  models  aimed  at  developing  a  suitable  theory 
of  the  zero-crossing  intervals  of  Gaussian  processes.  For 
example,  we  see  Immediately  that  a  theory  based  on  the 
Independence  of  successive  zero-crossing  intervals  is  destined 
to  be  inadequate. 
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Simllar  photographs  were  made  for  all  of  the  other  random 
processes  considered  in  this  report.  For  the  Gaussian  process 
having  .the  power  spectral  density  W^(f),  the  correlation 
coefficient  is  non-zero  when  n  s  4,  and  the  ith  and  (i+n  ',th  zero¬ 
crossing  intervals  are  statistically  dependent  when  n  s  4  .  For 
all  of  the  other  random  processes,  the  correlation  coefficient 
practically  vanishes  when  n  ^  3,  and  the  statistical  dependence 
between  the  ith  and  (i+n)th  zero-crossing  intervals  practically 
vanishes  when  n  2  4. 
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X.  APPLICATIONS 

The  sampling  technique  for  generating  Gaussian  noise  which 
was  described  in  Section  IIIB  is  one  application  resulting  from  this 
study  of  the  zero-crossing  intervals.  Another  application  concerns 
an  Integral  technique  for  measuring  pulse  duration  and  has  been 
reported  in  the  literature  (36). 

In  certain  countermeasure  applications  of  the  zero-crossing 
phenomenon  one  requires  a  modulating  noise  waveform  having 
an  "optimum"  Pq(t),  the  probability  density  function  for 
successive  zero-crossing  Intervals.  Usually  one  asks  that  a 
certain  mean,  Eq(t),  and  a  certain  standard  deviation,  be 
associated  with  the  Pq(t).  The  other  required  properties  of 
P^(t)  for  "optimum"  conditions  are  somewhat  vague.  Here  we 
shall  give  a  few  results  concerning  the  "optimum"  P^(t). 

A.  THE  MOST  RANDOM  DISTRIBUTION  OF 
SUCCESSIVE  ZERO-CROSSING  INTERVALS 
Consider  a  sequence  of  successive  zero-crossing 
intervale  {tj}  generated  by  a  sample  function  of  an  ergodic  random 
process  having  a  finite  expected  number  of  zeros  per  unit  time. 

Let  each  of  the  Tj  have  a  common  probability  density  Pq(t),  and 
lot  the  n  dimensional  joint  probability  density  for  the 
sequence  of  intervals  {tj}  be  denoted  by  P{j In 
order  to  determine  the  moot  random  distribution  of  successive 


I 
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zero-croBsing  intervals,  we  adopt  the  criterion  that  the  relative 
entropy  per  zero- crossing  interval,  H',  given  by 


,00  oo 

H' =  -  l^-.r  •••!  P(t^,T2.  .  .  .TjinP(T^,T2.  .  .  .T^)dT^.  .  .dr^ 


is  a  measure  of  the  randomness.  By  a  well-known 


(115) 


principle  from  information  theory,  the  maximum  of  H*  occurs 
when  the  sequence  of  successive  zero-crossing  intervals,  {t^},  are 
statistically  independent.  That  is  when 

P(ti,  T2.  •  • -Tn)  =  P^(t)  for  all  n  (116) 

Accordingly,  in  order  to  maximize  H'  we  need  only  msudmize 
the  relative  entropy,  H,  given  by 


(117) 


Let  the  constraints  be: 


J  P„(T)dT=l 


(118) 


tP„(t)  dr  =  E  (t) 


(119) 


o 


(120) 


By  applying  the  calculus  of  variations  we  find  that  the  maximum 
relative  entropy,  H,  results  when: 
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P 
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(t)  :: 


1 

Dv^  ®(^) 


(r-mf 

2 


where 


(121) 


E^(t)  =  m  +  XD 

E  (t^)  =  +  XDm  + 

o 


4(x)  = 


1 

■s/Zv 


2 

-X 


/2 


e 


»(x)  =  J  dy 

^2-11  -CD 
2 

Given  E^(t)  and  E^(t  )  one  can  determine  m  and  D  and  hence 
P^(t)  by  using  the  tables  reported  by  K.  Pearson  (37). 

Accordingly,  for  a  given  mean,  E^(t),  and  a  given  standard 

deviation,  ^  most  random  distribution 

of  Buccessuve  zero- crossing  intervals  is  represented  by  the 
truncated  Gaussian  probability  density  function  shown  in  Figure  67 
and  by  Equation  (116).  Since  successive  zero-crossing  intervals 
must  be  statistically  independent  in  order  to  yield  the  maximum 
relative  entropy  per  zero  crossing  interval,  the  theorem  of 
Section  IX  implies  that  a  Gaussian  process  or  an  infinitely  clipped 
Gaussian  process  can  not  generate  the  most  random  distribution 
of  successive  zero-crossing  intervals.  Moreover,  for  Gaussian 
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Figure  67  The  Truncated  Gaussian  Probability 
Density,  P  (t).  Represents  the  Most 
Random  Distribution  of  Successive 
Zero-Crossing  Intervals  Having  a 
Given  Mean  and  a  Given  Standard 
Deviation. 


T 
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processea  of  finite  memory,  Py{T)  for  large  r  must  have  an 
asymptotic  exponential  (5,  13,  38,  41)  behavior  such  as: 

PJt).  .  (122) 

Physically,  this  means  that  the  random  process  having  a  Pq(t) 
given  by  Equation  (121)  cannot  be  a  Gaussian  process  or  an 
infinitely  clipped  Gaussian  process.  That  is,  one  cannot  obtain 
a  random  process  having  a  Pq(t)  given  by  Elquation  (121)  by  the 
familiar  technique  of  spectrally  shaping  "white*  Gaussian  noise 
with  a  linear  filter. 

It  is  interesting  to  note  that  if  the  constraints  are  given  by 
Equations  (118)  and  (119)  only,  then  the  maximum  entropy  H 
occurs  when 

T 

The  Gaussian  process  having  power  spectral  density  Wj(f) 
generates  a  Pj,(t)  that  approximates  an  exponential  density  of  the  form 
given  by  Equation  (12  3)  as  is  seen  in  Figure  18.  Furthermore, 
successive  zero-crossing  intervals  of  this  Gaussian  process 
are  approximately  Independent  as  is  seen  in  Figure  19. 

B.  TECHNIQUE  FOR  ESTIMATING  THE  ASYMPTOTIC 

EXPONENTIAL  BEHAVIOR  OF  P^(t) 

Here  we  describe  a  technique  for  estimating  the  a 

in  Equation  (122). Lewie  (39)  and  McFadden  (32)  showed  that: 

E^(t^)  =  2E^(t)  J  ‘”p(O.T)dT  .  (124) 

o 


( 


-139- 

If  Pq(t)  ie  asymptotically  exponential  then  p(0,t)  and  Z(1,t)  are 
also  asymptotically  exponential  as  is  seen  from  Equations  39  and 
41.  Hence,  in  order  to  determine  a  in  Equation  (1Z2)  we  determine 
a  "reasonable"  exponential  to  tie  onto  the  measured  p(0,  t)  such 
that  Equation  (124)  is  satisfied  for  measured  E^(t^)  and  theoretical 
E^(t).  Semi-log  plots  of  Pq(t),  p(0,t)  and  Z(1,t)  are  useful  for 
determining  where  to  tie  on  the  exponential .  However,  the  point 
of  "tie  on"  is  still  somewhat  arbitrary.  Accordingly,  we  end  up 
with  at  most  rough  estimates  of  a.  Using  normalized  time,  some 
rough  estimates  of  associated  with  the  asymptotic  exponential 
behavior: 

are  given  below. 

For  tie  on  point  at  u^  =  8,  Uj  =  0.28 
For  tie  on  point  at  u^  =  12,  -  0.  1 1 

For  tie  on  point  at  u^  =  4,  a ^  =  0.47 
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XI.  CONCLUSIONS 

The  reaulta  of  thia  study  of  the  zero-croasing  intervala  of 
random  proceaaea  should  prove  useful  for  guiding  the  development 
of  some  new  techniques  which  are  useful  in  both  science  and 
technology. 

The  sampling  technique  for  generating  Gaussian  noise  is 
perhaps  the  simplest  technique  available  for  generating  such  noise. 

The  Integral  technique  for  measuring  pulse  duration  should 
be  especially  useful  for  airborne  or  space  instrumentation  applications. 

At  present  none  of  the  probabilities  or  probability  densities 
presented  in  this  report  can  be  derived  explicitly  by  analytical  methods. 

The  probability  density  Pq(t)  for  the  case  of  Gaussian  noise 

1  -t/tt 

having  a  power  spectral  density  Wj(f)  is  not  exactly  —  e  '  . 

The  hypothesis  that  successive  zero-crossing  intervals  of 
Gaussian  processes  from  a  wide  sense  Markov  chain  is  false  for 
some  Gaussian  processes. 

The  experimental  Initial  values  of  Pq(u^)  and  P|(u^)  for 
n  =  1,  3  obtained  by  extrapolation  agree  satisfactorily  with  the 
approximate  theoretical  initial  values. 

The  results  of  this  study  should  prove  useful  for  guiding  the 
development  of  models  aimed  at  developing  a  suitable  theory  for  the 
zero-crossing  Intervals  of  certain  random  processes.  In  general 
the  models  must  not  assume  that  successive  zero-crossing 
Intervals  of  random  processes  are  statistically  independent. 


For  the  Gauesian  process  having  the  power  spectral  density 
W^(f),  the  correlation  coefficient  is  nonzero  when  n  ^  4,  and 
the  ith  and  (i4-n)th  zero-crossing  intervals  are  statistically 
dependent  when  n  4.  For  all  of  the  other  random  processes, 
the  correlation  coefficient  practically  vanishes  when  n  a  3, 
and  the  statistical  dependence  between  the  ith  and  (i+n)th  zero¬ 
crossing  intervals  practically  vanishes  when  n  ^  4. 
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xn.  APPENDIX 

For  a  Gauatian  proceia  having  a  normalized  autocorrelation 
function  p(t).  the  conditional  probability,  U(t)  dr,  that  a  zero 
occura  between  t  +  t,  t  +  t  +  dr  given  a  zero  at  t  ia  given  by: 

U(T)dT  -  ^  [-p"(0)]  ^  -  p2(T)j  ^  l^l  +  Htan'^H  j 

(126) 

where  ^ - 

-  J  s  tan’^H  <  J 

2  ■  ■ 

M22*-p'’(0)  [l  -  p2(T)]  -[p'(t)] 

M2  3=p"(T)[i-pV)]  +  p(t)[p*(t)]^  . 

For  a  Gauasian  proceaa  having  a  normalized  autocorrelation 
function  p(t),  the  conditional  probability,  Q(t)  dr,  that  a  downward 
zero-croaaing  occura  between  t  +  t,  t  +  t  +  dr  given  an  upward 
zero-croaaing  at  time  t  ia  given  by: 

Q(T)dT=|l  1^.  p"(0)j’^  [^][^  -  P^(T)j’'^j^l  +  Hcot-^-H)j 

(127) 

where  — 

Os  cot  (-H)  S  IT 
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1 


The  normalized  autocorrelation  function,  r(u^,a),  for  the 
infinitely  clipped  process  ^)»  defined  by  Equation  (82),  is 
given  by: 


00  ,2 
h  ,  n 


I  -l^T 


n  =1 
odd  n 


p*^(u  )-8 
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oo  y  mu 
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oddm 


00  00  j^2 
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n=l  m=l 
odd  (m+n ) 


where 


m,  n 


(2m!) 


,n  m 
2  a 


2(!!2^;m+i:-al 

p2[2_-^m^)  J 


=  confluent  hypergeometric  function 
r  =  gamma  function 

a  =  signal-to-noise  power  ratio. 


The  normalized  autocorrelation  function,  r^(u^,  a),  for  the 
infinitely  clipped  process  t)(u,  a)t  defined  by  Equation  (106),  is 
given  by: 


r^(Uo.a)  = 
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where 
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